THE FINITE NON-PERIODIC TODA LATTICE: A GEOMETRIC AND 
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Abstract. In 1967, Japanese physicist Morikazu Toda published the seminal papers |78| and 
|79) , exhibiting soliton solutions to a chain of particles with nonlinear interactions between nearest 
neighbors. In the decades that followed, Toda's system of particles has been generalized in different 
directions, each with its own analytic, geometric, and topological characteristics that sets it apart 
from the others. These are known collectively as the Toda lattice. This survey describes and 
compares several versions of the finite non-periodic Toda lattice from the perspective of their 
geometry and topology. 
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1. Outline of the paper 
We organize the paper as follows: 

Section [2] introduces the real non-periodic Toda lattices. It begins with two formulations of the 
Toda lattice in which the flow obeys a Lax equation on a set of real tridiagonal matrices whose 
subdiagonal entries are positive. The matrices are symmetric in one formulation and Hessenberg in 
the other. The flows exist for all time and preserve the spectrum of the initial condition, and the 
system is completely integrable. The references cited for these results are [41 |2"T1 |2"51 1551 1551 1551 1551 175] , 
When the matrices of the Toda lattice system are extended to allow the subdiagonal entries to take 
on any real value, the two forms of the Toda lattice differ in the behavior of their solutions and 
in the topology of their isospectral sets. In the symmetric case, the flows exist for all time and 
the isospectral manifolds are compact, while in the Hessenberg form, the flows blow up in finite 
time and the isospectral manifolds arc not compact. For these results, the cited references are 
[T51 [T71 [501 15T1 [521 [551 [50] , Section 1 concludes with the full symmetric real Toda lattice, in which 
the flows evolve on the set of all real symmetric matrices. The symplectic structure comes from the 
Lie-Poisson structure on the dual of a Borel subalgebra of sl(n,M.). In the full Toda lattice, more 
constants of motion are needed to establish complete intcgrability. Cited references for this material 
are H2HHIIZ7]. 

Section [3J describes the iso-level sets of the constants of motion in complex Toda lattices. The 
complex tridiagonal Toda lattice in Hessenberg form differs from the real case in that the flows no 
longer preserve "signs" of the subdiagonal entries. The flows again blow up in finite (complex) time, 
and the isospectral manifolds are not compact. These manifolds are compactified by embedding 
them into a flag manifold in two different ways. In one compactification of isospectral sets with 
distinct eigenvalues, the flows enter lower-dimensional Bruhat cells at the blow-up times, where 
the singularity at a blow-up time is characterized by the Bruhat cell [25J [53]. Under a different 
embedding, which works for arbitrary spectrum, the Toda flows generate a group action on the flag 
manifold, where the group depends on how eigenvalues are repeated. The group is a product of a 
diagonal torus and a unipotent group, becoming the maximal diagonal torus when eigenvalues are 
distinct and a unipotent group when all eigenvalues coincide |70j . These group actions, together 
with the moment map [5 , 34 , 44 , are used in [72] to study the geometry of arbitrary isospectral sets 
of the complex tridiagonal Toda lattice in Hessenberg form. These compactified iso-level sets are 
generalizations of toric varieties |59j . Further properties of the actions of these torus and unipotent 
components of these groups are found in [301 ESI EH E2] [74j [75] . The survey then considers the full 
Kostant-Toda lattice, where the complex matrices in Hessenberg form are extended to have arbitrary 
complex entries below the diagonal. The techniques for finding additonal constants of motion in 
the real symmetric case are adapted to the full Kostant-Toda lattice to obtain a complete family 
of integrals in involution on the generic symplectic leaves; the geometry of a generic iso-level set 
is then explained using flag manifolds j26j [33] HQ] [67] . Nongeneric flows of the full Kostant-Toda 
lattice are described in terms of special faces and splittings of moment polytopes, and monodromy 
around nongeneric iso-level sets in special cases is determined [55J [55J [55J [70] . 

Section [4] provides other extensions of the Toda lattice. In this paper, we only discuss those related 
to finite non-periodic Toda lattices. We first introduce the flow in the Lax form on an arbitrary 
diagonalizable matrix, which can be integrated by the inverse scattering method (or equivalently 
by the factorization method) [52]. We then show how the tridiagonal Hessenberg and symmetric 
Toda lattices, which are defined on the Lie algebra of type A (that is, sl(n)), are extended to 
scmisimple Lie algebras using the Lie algebra splittings given by the Gauss decomposition and the 
QR decomposition, respectively []21 [I5J [T5J [5^ [35] [50] . As an important related hierarchy of flows, 
we explain the Kac-van Moerbeke system, which can be considered as a square root of the Toda 
lattice [35] [42] . We also show that the Pfaff lattice, which evolves on symplectic matrices, is related to 
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the indefinite Toda lattice [1. 48 , 49 . As another aspect of the Toda lattice, we mention its gradient 
structure [7J |5J EH HH H3] , which explains the sorting property in the asymptotic behavior of the 
solution and can be used to solve problems in combinatorian optimization and linear programming 

Section [5] explains the connections with the KP equation (4TJ [43] in the sense that the r-functions 
as the solutions of the Toda lattices also provide a certain class of solutions of the KP equation in a 
Wronskian determinant form [Bj |5TJ El EliSl 13 [37] . This is based on the Sato theory of the KP 
hierarchy, which states that the solution of the KP hierarchy is given by an orbit of the universal 
Grassmannian |64j . The main proposition in this section shows that there is a bijection between the 
set of r-functions that arise as a k x k Wronskian determiant and the Grassmannian Gr(fc,n), the 
set of all fc-dimensional subspaces of R". We present a method to obtain soliton solutions of the KP 
equation and give an elementary introduction to the Sato theory in a finite-dimensional setting. The 
moment polytopes of the fundamental representations of SL(n, C) play a crucial role in describing 
the geometry of the soliton solutions of the KP equation. 

Section [6] shows that the singular structure given by the blow-ups in the solutions of the Toda 
lattice contains information about the integral cohomology of real flag varieties [HI [19]. We begin 
with a detailed study of the solutions of the indefinite Toda lattice hierarchy [51] . The singular 
structure is determined by the set of zeros of the r-functions of the Toda lattice. First we note that 
the image of the moment map of the isospectral variety is a convex polytope whose vertices are the 
orbit of the Weyl group action [T51 [21] • Each vertex of the polytope corresponds to a fixed point of 
the Toda lattice, and it represents a unique cell of the flag variety. Each edge of the polytope can be 
considered as an orbit of the sl(2) Toda lattice (the smallest nontrivial lattice), and it represents a 
simple reflection of the Weyl group. There are two types of orbits, either regular (without blow-ups) 
or singular (with blow-ups). Then one can define a graph, where the vertices are the fixed points of 
the Toda lattice and where two fixed points are connected by an edge when the sl(2) flow between 
them is regular. If the flow is singular, then there is no edge between the two points. The graph 
defined in this way turns out to be the incidence graph that gives the integral cohomology of the 
real flag variety, where the incidence numbers associated to the edges are either 2 or —2 [45] [20], 
We also show that the total number of blow-ups in the flow of the Toda lattice is related to the 
polynomial associated with the rational cohomology of a certain compact subgroup [141 [T51 IT5] . 



2. Finite non-periodic real Toda lattice 

Consider n particles, each with mass 1, arranged along a line at positions g 1; q n . Between each 
pair of adjacent particles, there is a force whose magnitude depends exponentially on the distance 
between them. Letting p k denote the momentum of the kth particle, and noting that 4fQk = Pk 
since each mass is 1, the total energy of the system is the Hamiltonian 

1 n n — 1 

(2.1) H(p,q)= - £ pl+ £ e -(»+i-«») . 



2 



k=i fe=i 



The equations of motion 



(2.2) 



dg k _ dH_ 
dt dp k ' 

dpk _ dH 
dt dq k 
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(2.3) 



Pk 



yield the system of equations for the finite non-periodic Toda lattice, 

dqk 
dt 
dpk 
dt 

Here we set e~^ qi ~ q °^ = and e~^ qn+x ~ qn > = with the formal boundary conditions 



k = 1, n, 

e -(Sfc+i-<r*) + e -{lh-qk-i) k= i n 



q = -co, 



Qn+l 



In the 1970 's, the complete integrability of the Toda lattice was discovered by Henon and 
Flaschka |27j in the context of the periodic form of the lattice, where the boundary condition is 
taken as qo — q n . Henon [38 found analytical expressions for the constants of motion following 
indications by computer studies at the time that the Toda lattice should be completely integrable. 
That same year, Flaschka [27] [25J (independently also by Manakov [SB]) showed that the periodic 
Toda lattice equations can be written in Lax form through an appropriate change of variables. The 
complete integrability of the finite non-periodic Toda lattice was established by Moser [58] in 1980. 

A system in Lax form |55] gives the constants of motion as eigenvalues of a linear operator. In 
the finite non-periodic case, there are two standard Lax forms of the Toda equations. 



2.1. Symmetric form. Consider the change of variables (Flaschka [27J, Moser [55] ) 

1 

— < 
2 

1 



(2.4) 



_g-j(qfc+i-<?fc) 



:Pk, 



k 
k 



1 



In these variables, the Toda system (2.3 1 becomes 

da k 

(2.5) 



dt 
db k 
Itt 



a k (b k +i - b k ), 
2(« 2 fc -aLi), 



with boundary conditions 

This can be written in Lax form as 
(2.6) 
with 

(2.7) 



a = 0, a„ 



k = 1, n 
k = 1, n 

= . 



dt 



L(t) = [Skew(L(t)), L(t)] 



( 6 i 



V 



/ 



(2- 



Skew(L) := 



a i 



bn ) 



-ax 



-a„_i / 
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Any equation in the Lax form 4zL = [B,L] for matrices L and B has the immediate consequence 
that the flow preserves the spectrum of L. To check this, it suffices to show that the functions tr(L k ) 
are constant for all k. One shows first by induction that 4iL k = [B,L k ] and then observes that 



|[tr(L* 



tr[J^(L fe )] = tr[B, L k ] — 0. We now have n — 1 independent invariant functions 

1 



H k {L) = 



txL 



k+l 



The Hamiltonian (2.1) is related to Hi(L) by 

Hi(L) 



1 



with the change of variables (2.4 1 



If we now fix the value of H (L) = trL (thus fixing the momentum of the system), the resulting 
phase space has dimension 2(n — 1). With total momentum zero, this is 

' / h ai \ 



(2.9) 



S 



V 



>0, b t e 



o 



A property of real tridiagonal symmetric matrices (2.7 1 with ^ for all k is that the eigenvalues 
Afe are real and distinct. Let A be a set of n real distinct eigenvalues, and let M. = {L £ S : spec(i) = 
A }. Then S = Da-Ma- S is in fact a symplectic manifold. Each invariant function Hk(L) generates 
a Hamiltonian flow via the symplectic structure, and the flows are involutive with respect to the 
symplectic structure (see [1] for the general framework and [29 for the Toda lattice specifically) . We 



will describe the Lie-Poisson structure for the Toda lattice equations (2.6 1 in Section |2.5| however, 
we do not need the symplectic structure explicitly here. 

Moser [58j analyzes the dynamics of the Toda particles, showing that for any initial configuration, 
Qk+i ~ Qk tends to oo as t — * ±oo. Thus, the off-diagonal entries of L tend to zero as t -> ±oo so 
that L tends to a diagonal matrix whose diagonal entries arc the eigenvalues. We will order them 
as Ai < A2 < ••• < A n . The analysis in [58 shows that L(oo) = diag(A„, A n _i, • • • ,Ai) and 
L{— 00) = diag(Ai, A2, • • • , A„). The physical interpretation of this is that as t — > —00, the particles 
qk approach the velocities Pk(— 00) = — 2Afc, and as t — > 00, the velocities are interchanged so that 
Pfe(oo) = — 2A„_fc. Asymptotically, the trajectories behave as 



Pk(t) 



where A^" 



X k and A, 



A r , 



Symes solves the Toda lattice using matrix factorization, the QR-factorization; his solution, which 
he v erifies in |77j and proves in a more general context in |76j is equivalent to the following. To solve 
(2.5) with initial matrix L(0), take the exponential e tL ^ and use Gram-Schmidt orthonormalization 
to factor it as 



(2.10) 



JL(0) = fe( t ) r (f) 



where k(t) G SO(n) and r(t) is upper-triangular. Then the solution of (2.51 is 
(2.11) L(t) = fc- 1 (t)L(0)fc(t) = r(t) J L(0)r" 1 (i) . 

Since the Gram-Schmidt orthonormalization of e tL ^ can be done for all t, this shows that the 



solution of the Toda lattice equations (2.51 on the set (2.9 1 is defined for all t 
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We also mention the notion of the r-functions which play a key role of the theory of integrable 
systems (see for example |39j). Let us first introduce the following symmetric matrix called the 
moment matrix, 

(2.12) M(t) := e tL ^ T ■ e tL ^ = e 2tL ^ = r T (t)k T (t)k(t)r(t) = r T (t)r(t) , 

where r T denotes the transpose of r, and note k T = k^ 1 . The decomposition of a symmetric matrix 
to an upper-triangular matrix times its transpose on the left is called the Cholesky factorization. 
This factorization is used to find the matrix r, and then the matrix k can be found by k — e^C)^ 1 . 
The T-functions, Tj for j — 1, . . . ,n — 1, are defined by 



(2.13) 



r 3 (t) := det (M j (t)) = l[r i (t) 2 , 



where Mj is the j x j upper- left submatrix of M, and we denote diag(r) = diag(r x . . . , r n ). Also 
note from (2.111, i.e. L(t)r(t) — r(t)L(0), that we have 

r k +i(t) 



aj (t) = aj{0) 



With ( |2.13[ ), we obtain 
(2.14) 



GjO) = a j(°) 



r k (t) ' 
v /r 3+1 (^)r 3 _ 1 (t) 

Tj(t) 



and from this we can also find the formulae for bj(t) of L(t) as 



2dt 



In 



Tj-l(t) 



One should note that the r-functions are just defined from the moment matrix M = e 2tLly0 \ and 
the solutions (aj{t),bj(t)) are explicitly given by those r-functions without the factorization. 



2.2. Hessenberg form. The symmetric matrix L in (2.7), when conjugated by the diagonal matrix 
, a„_i), yields a matrix Y = DLD^ 1 in Hessenberg form: 



D = diag(l,ai, 
(2.15) 



Y 



( h 1 



V 



\ 



b n j 



The Toda equations now take the Lax form 
(2.16) 
with 



dt 



2[Y,F) 



F = 



/ 



V 



o 
o / 



Notice that Y{t) satisfies ( |2.16[ ) if and only if X(t) = 2Y(t) satisfies 

d 



(2.17) 



dt 



X=[X,U N _X) 
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where H_\f_A is the strictly lower-triangular part of A obtained by setting all entries on and above 



the diagonal equal to zero. Equation (2.171 with 



/ h 



(2.18) 



X 



1 



91 



V 



9n 



fn j 



is called the asymmetric, or Hessenberg, form of the non-periodic Toda lattice. Again, since the 
equations are in Lax form, the functions Hf~(X) = jrj-j- trX k+1 . are constant in t. 



Notice that the Hessenberg and symmetric Lax formulations of (2.3 1 are simply different ways 
of expressing the same system. The solutions exist for all time and exibit the same behavior as 
t — > ±oo in both cases. However, when we generalize the Toda lattice to allow the subdiagonal 
entries to take on any real value, the symmetric and Hessenberg forms differ in their geometry and 
topology and in the character of their solutions. 

2.3. Extended real tridiagonal symmetric form. Consider again the Lax equation 
(2.19) jL= [Skew(L), L] , 



where we now extend the set of matrices L of the form (2.7) by allowing a*; to take on any real 
value. (Recall that in our original definition of the matrix L, each ak was an exponential and was 
therefore strictly positive.) As before, bk may be any real value, and SkewL is defined as in (2.8 1. 



Given a ny i nitial matrix in this extended form, the factorization method of Symes, described 



in Section 



2.1 



still works. Indeed, for any such initial matrix £(0), e tL ' ' can be factored into 
(orthogonal) x (upper-triangular) via the Gram-Schmidt procedure, i.e. the QR-factorization, and 
one can verify as before that L(t) — k~ 1 (t)L(0)k(t), where k(t) is the orthogonal factor. Thus, in 
the extended symmetric form, solutions are still defined for all t. Given the initial eigenmatrix of L, 



an explicit solution of (2.191 can be found in terms of the eigenmatrix of L(0) using the method of 
inverse scattering. 

For a general Lax equation ^Y — [B,Y], if Y(0) has distinct eigenvalues Ai,...,A„, then the 
inverse scattering scheme works as follows. Let A be the diagonal matrix of eigenvalues, A = 
diag(Ai, A„), and let $(t) be a matrix of normalized eigenvectors, varying smoothly in t, where 
the fcth column is a normalized eigenvector of Y(t) with eigenvalue A&. Then the Lax equation 
4tY = [B,Y] is the compatibility condition for the equations 



(2.20) 
(2.21) 

We see this as follows: Denoting (•)' 
(y$)' = ($A)' 



Y(tWt) 
d_ 
dt 



*(*) = B(t)*(t) 



(F)'$ + y($)' = ($)'a 

(Y)'<S> + YB(j) = B$A 
(Y)'$= -YB<$> + BY<$> 
(y)'$ = [B,Y]<f> 
(Y)' = [B,Y] . 



by (2.211 



by (23U| 



The inverse scattering method solves the system (2.201 and (2.21) for <p(t) and then recovers Y(t) 



from (2.21 1. Since B(t) is defined as a projection of Y(t), which can be written in terms of and 
A, we obtain a differential equation for $(t) by replacing B(t) in (2.211 by its expression in terms 
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(1,2,3) 

(2,1 ,3) jf^^Sc > 3 ' 2 ) 




(2,3,1) "S^v^ (3,1,2) 
(3,2,1) 



(M + . 



(M. 



(M. 



FIGURE 2.1. The Tomei manifold .Ma for the sZ (3,1ft) Toda lattice. The 3-tuples 
(i,j,k) on the vertices indicate the diagonal matrices L = diag(A^, X -,, Afc). Each 
hexagon M €l ,e 2 corresponds to the moment polytope (see Section 3.1 I for the Toda 
lattice with the signs (ex, e^) = (sgn(ai), sgn(a2)). The boundaries correspond to 
the sZ(2,R) Toda lattices associated with either a\ = or a-i — 0. The Tomei 
manifold A4\ is given by gluing edges of the hexagons. For example, the edge 

between (1,2,3) and (1,3,2) in M++ is glued with the same edge in M h since 

this edge indicates ai = 0. The other gluing shown in the figure is for the edges 
corresponding to 02 =0. 



of <&(i) and A. Given Y(0), we then obtain $(0) from (2.20) with t = 0. This solution is in fact 
equivalent to the QR factorization given above (see Section 2.1|. 

For real matrices of the form (2.7), the inverse scattering method can be used on the open dense 
subset where all are nonzero. This is because a real matrix L of the form (2.7 1 has distinct real 
eigenvalues if 7^ for all k. The eigenvalues are real because L is a real symmetric matrix; the 
fact that they are distinct follows from the tridiagonal form with nonzero a^, which forces there to 
be one eigenvector (up to a scalar) for each eigenvalue. 

Let M\ denote the set of n x n matrices of the form (2.7 1 with fixed eigenvalues Ai < A2 < 
• • • < A n . .Ma contains 2" _1 components of dimension n — 1, where each component consists of all 
matrices in .Ma with a fixed choice of sign for each a^. The solution of (2.19) with initial condition 
in a given component remains in that component for all t; that is, the solutions preserve the sign 
of each et^. Each lower-dimensional component, where one or more is zero and the signs of the 
other a,k are fixed, is also preserved by the Toda flow. Adding those lower dimensional components 
gives a compactification of each component of .Ma with fixed signs in a^'s. Tomei [5U] shows that 
A4\ is a compact smooth manifold of dimension n — 1. In the proof of this, he uses the Toda flow 
to construct coordinate charts around the fixed points. Tomei shows that Ma is orientable with 
universal covering IR™ -1 and calculates its Euler characteristic. 

In his analysis, Tomei shows that Ma contains 2 n ~ 1 open components diffcomorphic to K" _1 . On 
each of these components, ^ for all k, and the sign of each au is fixed. They are glued together 
along the lower-dimensional sets where one or more is zero. For example, in the case n = 3, there 

are four 2-dimensional components, denoted as M++ 7 M^ ,M |_, and M , according to the signs 

of ai and ai. The closure of each component is obtained by adding six diagonal matrices where all 
the cik vanish (these are the fixed points of the Toda flow) and six 1-dimensional sets where exactly 
one dfc is zero. Denote the closure of M++ by M++, and so on. The boundary of M++, for example, 
contains three 1-dimcnsional sets with a\ = and 02 > 0. Each is characterized by having a fixed 
eigenvalue as its first diagonal entry. The other three 1-dimensional components in the boundary of 
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A4++, have a\ > and a 2 — 0, with a fixed eigenvalue in the third diagonal entry. Adding those 
boundaries with the 6 vertices corresponding to the diagonal matrices L — diag(Ai, A,-, A&) gives the 
compactified set A4++. Notice that the three 1-dimcnsional components with a\ — 0,a 2 > 0, and a 

fixed Afc in the first diagonal entry also lie along the boundary of A4 |_; the other three components, 

with di > and 02 = are shared by the boundary of AA++. In this manner, the four principal 
components are glued together along the subset of M\ where one or more vanish. In Figure 



2.1 we illustrate the compactification of the Tomei manifold A^a for the s/(3,IR) symmetric Toda 
lattice, 

M A = M ++ U M + _ U M- + U M— , 
where the cups include the specific gluing according to the signs of the as explained above. The 
resulting compactified manifold At a is a connected sum of two tori, the compact Riemann surface 
of genus two. This can be easily seen from Figure |2~T| as follows: Gluing those four hexagons, A4a 
consists of 6 vertices, 12 edges and 4 faces. Hence the Euler characteristic is given by x(-^a) = 
6 — 12 + 4 = —2, which implies that the manifold has genus g — 2 (recall x — 2 — 2g). It is also 
easy to see that A1a is orientable (this can be shown by giving an orientation for each hexagon so 
that the directions of two edges in the gluing cancel each other). Since the compact two dimensional 
surfaces are completely characterized by their oricntability and the Euler characters, we conclude 
that the manifold A^a is a connected sum of two tori, i.e. g = 2. 

The Euler characteristic of A4a (for general n) is determined in [SUj as follows. Let L = 
diag(A CT ( 1 ), A cr („)) be a diagonal matrix in A4\, where a is a permutation of the numbers {1, n}, 
and let r(L) be the number of times that a(k) is less than a(k + 1). Denote by E(n, k) the number 
of diagonal matrices in Ma with r(L) = k. Then the Euler characteristic of Ma is the alternating 
sum of the E(n, k): 

X(M K )= ±(-\) k E(n,k) . 

k=0 

[In M. Davis et al extends Tomei's result....] 

If the eigenvalues of the tridiagonal real matrix L are not distinct, then one or more must be 
zero. The set of such matrices with fixed spectrum where the eigenvalues are not distinct is not a man- 
ifold. For example, when n = 3 and the spectrum is (1, 1,3), the isospectral set is one-dimensional 
since one is zero. It contains three diagonal matrices Di — diag(3, 1, 1), D 2 — diag(l,3, 1), and 
D3 = diag(l, 1,3), and four 1-dimensional components. Each 1-dimensional component has a 1 in 
cither the first or last diagonal entry and a 2 x 2 block on the diagonal with eigenvalues 1 and 3, 
where the off-diagonal entry is either positive or negative. The two components with the lxl block 
in the last diagonal entry connect D\ and D2, and the two components with the lxl block in 
the first diagonal entry connect D 2 and D 3 . In Figure |2.2[ we illustrate the isospectral set of those 
matrices which is singular with a shape of figure eight. 

2.4. Extended real tridiagonal Hessenberg form. We now return to the Hessenberg form of 
the Toda equations with 

/ h 1 



(2.22) X = 



.91 



\ ffn-1 fn / 



as in (2.181, and allow the g/. to take on arbitrary real values. On the set of tridiagonal Hessenberg 



matrices X with and gu real, the Toda flow is defined by 
(2.23) j t X{t) = [X{t)^ N _X{t)} 
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(1,1,3) 



(3,1,1) 



(1,1,3) 



(1,1,3) (1,1,3) 

(0+) \(o-) 



< 1 - 3 ' 1 >\(+0) PAD \(+0) 0,3,1) \(_o) d,3,1) \(_ 0) 



\^+0) PA1) \^+ 



(3,1,1) 



(3,1,1) 



(3,1,1) 



(M 4 



(M. 



Figure 2.2. The isospectral set of the matrices having the eigenvalues {1,1,3}. 
Each polytope M. tl ^ t2 contains the set of matrices with either a\ = or a 2 = and 
the signs = sgn(a;). These polytopes are obtained by squeezing the polytopes for 



the semisimple sl(3) Toda lattice in Figure 2.1 according to the degeneration of the 
eigenvalues 2 — > 1. The gluing pattern is the same as in Figure 2.1 that is, identify, 

for example, the edge (0+) in with the same one in A4 h . The resulting 

variety is singular and has a shape of figure eight. 



as in ( |2.17 l. 

Recall that in the formulation of the original Toda equations, all the gk were positive, so that the 
eigenvalues were real and distinct. When for some k, the eigenvalues may now be complex 
or may coincide. We will see that this causes blow-ups in the flows so that the topology of the 
isospectral manifolds are very different from the topology of the Tomei manifolds described in the 
previous section. 

The matrices of the form (2.221 with gk ^ for all k are partitioned into 2™" 1 different Hamil- 
tonian systems, each determined by a choice of signs of the gk- Letting sj, = ±1 for k = 1, n and 
taking the sign of gk to be SkSk+i, Kodama and Ye |51j give the Hamiltonian for the system with 
this choice of signs as 



(2.24) 



where 



fk 



1 11 
H=IZ yl 

z k=l 



Vk 



SfeSfe+ie 

fc=i 



-(Xk + 1-Xk) 



k = 1, 



9k 



■ SkSk+i e 



k = 1, 



1 



The system (2.231 is then called the indefinite Toda lattice. The negative signs in (2.241 correspond 
to attractive forces between adjacent particles, which causes the system to become undefined at 
finite values of t, as is seen in the solutions obtained by Kodama and Ye in [50] and [51] by inverse 
scattering. 

The blow-up s in t he solutions arc also apparent in the factorization solution of the Hessenberg 
form. To solve (2.23) with initial condition Y(0), factor the exponential e tx ^ as 

(2.25) e* x(0) = n(t)b(t) , 

where n(t) is lower unipotent and b(t) is upper-triangular. Then, as shown by 



and [5T] . 



(2.26) 



X(t) = n" 1 (t)X(0)n(t) = b^X^b-^t) 
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solves (2.23). Notice that the factorization (2.251 is obtained by Gaussian elimination, which multi- 
plies e txi ^^ on the left by elementary row operations to put it in upper-triangular form. This process 
works only when all principal determinants (the determinants of upper left k x k blocks, which are 
the r-functions defined below) are nonzero. At particular values of t e R, this factorization can fail, 
and the solution (2.26) becomes undefined. 

The solutions (/&, <?fc) can be expressed in terms of the r-functions which are defined by 



(2.27) 



r k (t) :=det [(e tx (°)) fc ] =*[[<!&), 



where (e^'- -')^ is the k x k upper-left submatrix of e tL ^°\ and diag(6) = diag(<ii, . . . , d n ). With 
(2.261, we have 



(2.28) g k (t) = g k (0) 
The function fk{i) are given by 

(2.29) f k {t) 



dk+i(t) 
dk{t) 



9k(P) 



Tfe- 



-lft)Tfc-l(t) 



dt 



In 



Tfe(t) 
Tk-l(t) 



Now it it clear that the factorization (2.25 1 fails if and only if Tfc(t) = for some k. Then a blow-up 
(singularity) of the system (2.23 1 can be characterized by the zero sets of the r-functions. 

Example 2.1. To see how blow-ups occur in the factorization solution, consider the initial matrix 



When t ± -1, 



JX 



1+t 
-t 



t 

1-t 



1 

-t 
1+t 



1 + t 





t 

1 

1+t 



and the solution evolves as in (2.261. The r-function is given by T\{t) = 1 + t, and when t = — 1, 
this factorization does not work. However, we can multiply e~ ix " on the left by a lower unipotent 
matrix nT 1 (in this case the identity) to put it in the form wb, where w is a permutation matrix: 



-ix 



1 



-1 




This example will be taken up again in Section |3.2| where it is shown how the factorization using a 
permutation matrix leads to a compactification of the flows. 



In general, when the factorization (2.251 is not possible at time t — t, e tx ^ can be factored 
as e*^ -* = n(t) w b(t), where w is a permutation matrix. Ercolani, Flaschka, and Haine use 
this factorization to complete the flows (2.261 through the blow-up times by embedding them into 
a flag manifold. The details will be discussed in the next section, where we consider the complex 
tridiagonal Hessenberg form of the Toda lattice. 

Kodama and Ye find explicit solutions of the indefinite Toda lattices by inverse scattering. Their 
method is used to solve a generalization of the full symmetric Toda lattice in [SU] and is specialized 
to the indefinite tridiagonal Hessenberg Toda lattice in j 
and Ye make the change of variables 



For the Hamiltonian (2.24), Kodama 



1 



-(x k + 1 ~x k )/2 



(2.30) 



fc = l,. 



1. 
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together with t — ► It so that Hamilton's equations take the form 



(2.31) 



da k 
~dT 
db k 
dt 



2 a fe( s fe+i^fe+i — s kbk) 



Sk+ial - Sfe-iafe.! 



with ao = a n — 0. Here we switched the notation and b k from the original one in |50l 151] . This 
system is equivalent to (2.231 with fk = Skbk and g k = SkSk+iaf,. The system (2.311 can then be 
written in Lax form as 



(2.32) 



±L=[B,L], 



where L is the real tridiagonal matrix 

/ si&i s 2 a\ 

S\a\ S 2 &2 S 3 a 2 

(2.33) L = 





V o 



§71— lbji — 1 Sn^n—l 



and B is the projection 
(2.34) 



B = -[(L) >0 -(L) <0 ] 



The inverse scattering scheme for (2.34) is 
(2.35) 



dt 



$ = 5$ 



where A = diag(Ai, A„) and $ is the eigenmatrix of L, normalized so that 
(2.36) $£ -1 $ T = £ -1 , $ T S* = S 



with 5 = diag(si, s n ). Note that the matrix L of (2.33 1 is expressed as L = LS with the symmetric 



matrix L given by ( |2.7[ ) for the original Toda lattice. When 5 is the identity (2.361 implies that $ 
is orthogonal, and for S = diag(l, 1, — 1, —1), $ is a pseudo-orthogonal matrix in 0(p, q) with 
p + q = n. This then defines an inner product for functions / and jona set C, 

n 

(f,g) :=ff(\)g(X)d(x(X) = £ f(X k )g(X k )s^ , 

C k=l 
n 

with the indefinite metric d[i := £ s / 7 1 (5(A — Afc) d\ and the eigenvalues A& of L. Then the entries 

k=l 

of L can be expressed in terms of the eigenvector 0(A) = (</>i(A), . . . , <fi n (\)) T , i.e. L<fi(A) = X(p(X), 



(2.37) 



fe=i 



The explicit time evolution of $ can be obtained using an orthonormalization procedure on functions 
of the eigenvectors that generalizes the method used in [47^ to solve the full symmetric Toda hierarchy. 
A brief summery of the procedure is as follows: First consider the factorization (called the HR- 
factorization) , 



(2.38) 
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where r(t) is a lower triangular matrix and h(t) satisfies h T Sh = S (if S = I, then h 6 SO(n), 
i.e. the factorization is the QR-type). Then the eigenmatrix <fr(t) = (<f>i(\j))i<i j< n is given by 
$(t) = /i(t)$(0). Now one can write $(f) as 

$(t) = /i(i)*(0) = r _1 (t)r(t)fe(t)*(0) 

= r- 1 (t)e* x(0) $(0) = r- 1 (t)$(0)e* A . 

Since r(i) is lower triangular, this implies 

<f>%(\ t) = S P an R {0°(A)e At , . . . , 0°(A)e At } , i = 1, . . . ,n. 

Then using the Gram-Schmidt orthogonalization method, the functions <j>i(\, t) can be found as |51j . 



(2.39) 



<j>i{\t) 



-M 



where #(A) = ^(A,0), c l3 {t) 



SiCu ••• 

S1C 2 1 ••• 

SlCa • • • 
and Dk(t) 



-1C1,., 
-ic 2 ,., 



</>?(A) 



St-iCi,i_i $(A) 

|(SiCijft))i<i,j< fc |. 



inverse scattering problem (2.351 is then obtained from (2.391 using (2.371. 



The solution of the 
The matrix M := 



{cij)i<i,j<n is the moment matrix for the indefinite Toda lattice which is defined in the similar way 
as ( |2.12| ) Ti e. 

M(t) := e i i ^S- 1 ei L '^ T = ^e^^ 1 ^ , 
where we have used L(0)<I>o = and ^S^o = S. Then the r-functions are defined by 

(2.40) f k (t) = det (M fc (t)) = |(cii(t))i<«<*| = D k (t) . 

Sf ■■ Sk 



From (2.391 it follows that when fk(t) = for some k and time i, L(t) blows up to infinity as 
t — * t. In |51j . Kodama and Ye characterize the blow-ups with the zeros of r-functions and study the 
topology of a generic isospectral set .Ma of the extended real tridiagonal Toda lattice in Hessenberg 
form. 

It is first shown, using the Toda flows, that because of the blow-ups in L, A4\ is a noncompact 
manifold of dimension n — 1. The manifold is then compactified by completing the flows through 
the blow-up times. The 2x2 case is basic to the compactification for general n. The set of 2 x 2 
matrices with fixed eigenvalues Ai and A2, 



(2.41) 



Mi 




consists of two components, M+ with gi > and M- 

Ai 1 



A 2 



and L 



■ Ai < A 2 I 

with g\ < 0, together with two fixed points, 

A 2 1 
Ai 



Writing / 2 = Ai+A 2 — fi and substituting this into the equation for the determinant, /i/ 2 — <?i = AiA 2 , 
shows that M\ is the parabola 

(2.42) ff l = -(/i-Ai)(/i-A a ). 

This parabola opens down, crossing the axis g\ — at fi — Ai and fi — A 2 , corresponding to the 
fixed points L\ and L 2 . For an initial condition with g\ > 0, the solution is defined for all t; it flows 
away from p 2 toward p\. This illustrates what is known as the sorting property, which says that 
as t — ► 00, the flow tends toward the fixed point with the eigenvalues in decreasing order along the 
diagonal. The component with gi < is separated into disjoint parts, one with fi < Ai and the 
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(1,2,3) 

(2,1,3) (1,3,2) 




(2,3,1) ^ ^ (3,1,2) 
(3,2,1) 




(M. 



Figure 2.3. The compactification of the isospectral manifold A4\ for the indefinite 
sl(3, M.) Toda lattice. As in Figure 2.1 each hexagon indicates the moment polytope 
associated with the indefinite Toda lattice. The signs (ei, £2) in M. tl ^ are those of 
((71,(72) for t — ► —00, and each sign in the hexagons indicates the signs of (31,32)- 
The gluing rule according to the sign changes of gt is the same as that in the Tomci 
manifold, but the pattern is now different. For example, the edge between (2,3,1) 

and (3,2, 1) in M is now glued with that in M.^ The solid and dashed lines 

in the hexagons show the points where the solutions blow-up, i.e. = for k = 1 
(solid) and for k — 2 (dashed). Also the numbers in the hexagons indicate the 
number of blow-ups along the flow from t = —00 to +00 (see Section 6.3 1. 



other with fi > A 2 . The solution starting at an initial matrix with fi > A 2 flows toward the fixed 
point L2 as t — *■ 00. For an initial matrix with fi < Aj., the solution flows away from L±, blowing up 
at a finite value of t. By adding a point at infinity to connect these two branches of the parabola, 
the flow is completed through the blow-up time and the resulting manifold is the circle, S 1 . 

For general n, the manifold Ma with spectrum A contains n\ fixed points of the flow, where 
the eigenvalues are arranged along the diagonal. These vertices are connected to each other by 
incoming and outgoing edges analogous to the flows connecting the two vertices in the case n = 2. 
On each edge there is one gt that is not zero. As in the case n = 2, edges in which the blow-ups 
occur are compactified by adding a point at infinity. Kodama and Ye then show how to glue on 
the higher-dimensional components where more than one g^ is nonzero and compactify the flows 
through the blow-ups to produce a compact rt-dimensional manifold. The result is nonoricntable for 
n > 2. In the case n = 3, it is a connected sum of two Klein bottles. In Figure 2.3 we illustrate the 
compactification of Ma for the sZ(3,K) indefinite Toda lattice. With the gluing, the compactified 
manifold A^a has the Euler characteristic x(A^a) = —2 as in the case of the Tomei manifold (see 
Figure 2.1 1. The non-orientability can be shown by non-cancellation of the given orientations of the 
hexagons with this gluing. 

The compactification was further studied by Casian and Kodama [TS] (also see [T7]), where 
they show that the compactified isospectral manifold is identified as a connected completion of the 
disconnected Cartan subgroup of G = Ad(SL(n,M) ). The manifold is diffeomorphic to a toric 
variety in the flag manifold associated with G. They also give a cellular decomposition of the 
compactified manifold for computing the homology of the manifold. We will show more details in 
Section 6.3 where the gluing rules are given by the Weyl group action on the signs of the entries gj 
of the matrix X. 
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2.5. Full Symmetric real Toda lattice. We now return to the symmetric Toda equation 

(2.43) — L = [Skew(L), L] 

at 



as in (2.6 1, where L is now a full symmetric matrix with distinct eigenvalues. As in the tridiagonal 
case, Skew(L) is the skew-symmetric summand in the decomposition of L into skew-symmetric plus 
lower-triangular. Deift, Li, Nanda, and Tomei [22] show that (2.431 remains completely integrable 
even in this case. They present a sufficient number of constants of motion in involution and construct 
the associated angle variables. 

The phase space of the full symmetric real Toda lattice is the set of symmetric matrices, which 
we denote 

Sym(n) := {L e sl(n,R) : L T = L) . 

One can then define the Lie-Poisson structure (the Kostant-Killilov 2-form) on this phase space as 
follows. First we define a nondegenerate inner product (the Killing form), (A,B) = tr(AB), which 
identifies the dual space s/*(n,K) with sl(n,R). Then we consider the Lie algebra splitting, 

sl(n,R) = B- so(n) = B*_ ® so*{n) , 

where B_ is the set of lower triangular matrices. With the inner product, we identify 

B*_ S so X (n) = Sym(n), so*(n) = B± = JV_ , 

where A/"_ is the set of strictly lower triangular matrices, and A indicates the orthogonal comple- 
ment of A with respect to the inner product. The Lie-Poisson structure is then defined as follows: 
For any functions / and h on B*_ = Sym(n) , define 

{f,h}(L) = (L,[n B _Vf,u B _vh]) , 

where (X, V/) = lim e ^o £f(L + £ A), and Hjs_X is the projection of X onto B-. The Toda lattice 
(2.431 can now be expressed in Hamiltonian form as 

jh = {Hx,L}(L) with H X {L) = ^tr(L 2 ) . 



Using the Poisson structure, we can now extend equation (2.431 to define the Toda lattice hierarchy 
generated by the Hamiltonians Hk(L) = -j^tv(L k+1 ): 

(2.44) ^-L = {H k ,L}(L) = [Skew(L k ),L] with H k = ^tr(L fe+1 ) , 

OTfe K + 1 

where Skew(L fe ) = n so ( n ) VHk- The flow stays on a co-adjoint orbit in Sym(rt) = B*_. The Lie- 
Poisson structure is nondegenerate when restricted to the co-adjoint orbit, and the level sets of the 
integrals found in |22j are the generic co-adjoint orbits. 

Deift and colleagues find the constants of motion by taking the matrices (L)k obtained by removing 
the first k rows and last k columns of L. A co-adjoint orbit is obtained by fixing the trace of each 
(L)k- The remaining coefficients of the characteristic polynomials of (L)j. for < k < [n/2] (that is, 
all coefficients except for the traces) provide a family of [n 2 /4] constants of motion in involution on 
the orbit. Generically, (L)& has n — 2k distinct eigenvalues Ai^, \ n ~2k,k- The constants of motion 
may be taken equivalently as the eigenvalues A r .fc for < k < [\{n — 1)] and 1 < r < n — 2k. In this 
case the associated angle variables are essentially the last components of the suitably normalized 
eigenvectors of the 

In [27], Kodama and McLaughlin give the explicit solution of the Toda lattice hierarchy (2.441 
on full symmetric matrices with distinct eigenvalues by solving the inverse scattering problem of the 



16 



YUJI KODAMA 1 AND BARBARA SHIPMAN 2 



system 

L<S> = $A, 

= Skew(L fc )$ . 

with A = diag(Ai, . . . , A„). Since L is symmetric, the matrix $ of eigenvectors is taken to be 
orthogonal: 

L = $ A $ T 

with $ = [</>(Ai), </>(A n )], where the 4>(\k) 1S the normalized eigenvector of L wit h eige nvalue A&. 



The indefinite extension of the full symmetric Toda lattice (where L — LS as in (2.331 is studied 
in [37], where explicit solutions of <f>(\k,i) are obtained by inverse scattering. The authors also give 
an alternative derivation of the solution using the factorization method of Symes [77], where e tL ^ is 



factored into a product of a pseudo-orthogonal matrix times an upper triangular matrix as in (2.381 
(the HR- factorization) . 

3. Complex Toda lattices 

Here we consider the iso-spectral varieties of the complex Toda lattices. In order to describe 
the geometry of the iso-spectral variety, we first give a summary of the moment map on the flag 
manifold. The general description of the moment map discussed here can be found in |441 134] . 



3.1. The moment map. Let G be a complex semisimplc Lie group, H a Cartan subgroup of G, and 
B a Borel subgroup containing H. If P is a parabolic subgroup of G that contains J5, then G/P can 
be realized as the orbit of G through the projectivized highest weight vector in the projectivization, 
P(V), of an irreducible representation V of G. Let A be the set of weights of V, counted with 
multiplicity; the weights belong to TC^., the real part of the dual of the Lie algebra TC of H . Let 
{v a : a E A} be a basis of V consisting of weight vectors. A point [X] in G/P, represented by X e V, 
has homogeneous coordinates n a (X), where X — E Tr a (X)v a . The moment map as defined in |44] 

sends G/P into H^: 

A 4 : G/P ► H* R 



(3.1) 

[X] 



£ Wa(X)\ 2 <* 

a<£A 



Its image is the weight polytope of V, also referred to as the moment polytope of G/P. 

The fixed points of H in G/P are the points in the orbit of the Weyl group W through the 
projectivized highest weight vector of V; they correspond to the vertices of the polytope under the 
moment map. Let H ■ [X] be the closure of the orbit of H through [X] . Its image under fi is the 
convex hull of the vertices corresponding to the fixed points contained in H ■ [Xj; these vertices are 
the weig hts {a € W ■ a v : tt q (X) ^ 0}, where a v is the highest weight of V |5j. In particular, the 
image of a generic orbit, where no 7r Q vanishes, is the full polytope. The real dimension of the image 
is equal to the complex dimension of the orbit. 

In the case that P = B. V is the representation whose highest weight is the sum of the fundamental 
weights of G, which we denote as 5. Let v$ be a weight vector with weight S. Then the action of 
G through [vg] in P(V) has stabilizer B so that the orbit G • [vg] is identified with the flag manifold 
G/B. The projectivized weight vectors that belong to G/B are those in the orbit of the Weyl group, 
W = N(H)/H, through [vg], where N(H) is the normalize!' of H in G. These are the fixed points 
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FIGURE 3.1. The moment polytope of G/B for G = SX(4,C). The number set 
*i*2*3*4 indicates the weight L = i\L\ + ^2^2 + *3^3 + «4i4- The highest weight 
corresponds to the set 3210. 



of H in G/B. The stabilizer in W of [vg] is trivial so that in G/B, the fixed points of H are in 
bijection with the elements of the Weyl group. 

Now take G = SL(n, C), B the upper triangular subgroup, and H the diagonal torus. The choice 
of B determines a splitting of the root system into positive and negative roots and a system A of 
simple roots. The simple roots are Li — where i = 1, ...,n — 1 and L t is the linear function 

on 7i that gives the zth diagonal entry; the Weyl group is the permutation group £„, which acts by 
permuting the Li. 7ijj is the quotient of the real span of the Li by the relation L\ + • • • + L n = 0. 
7ijj may be viewed as the hyperplane in K™ where the sum of the coefficients of the Li is equal to 
1 + 2 + • • • + (n — 1). Let ...,«„} = {0, ...,n}. The moment polytope is the convex hull of the 
weights L = i 1 L 1 + i 2 -^2 + ''' +'ninj where (ii,«2j • • • i *n-ii *'«.) — (n — 1, n — 2, • • ■ ,1,0) corresponds 
to the highest weight. In Figure |3~Tj we illustrate the moment polytope for the flag manifold G/B 
of G = SL(4,C)- 

Let S be the set of reflections of Tt'^ in the hyperplanes perpendicular, with respect to the Killing 
form, to the simple roots (these are the simple reflections) . The group of motions of 7i Jj generated 
by S is isomorphic to W\ it is also denoted as W and referred to as the Weyl group of G. The 
vertices of the moment polytope are the orbit of W through 8. For w £ W, the moment map /i 
sends w[vs] to the vertex w(S). An arbitrary w G W can be written as a composition w = s r . . . s\ 
of simple reflections The length, l(w), of w with respect to the simple system A is the smallest 
r for which such an expression exists. 



3.2. Complex tridiagonal Hessenberg form. 

Hessenberg matrices 

/ A 1 



Here we consider the set M. of complex tridiagonal 



(3.2) 



X = 



.91 



V 



1 

fn J 



where the fk and gk are allowed to be arbitrary complex numbers. As before, the Toda flow is 



defined by (2.231 and the eigenvalues (equivalently, the traces of the powers of X) are constants of 
motion. The Hamiltonian 
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generates the flow 
(3.3) 



ax 

dt k 



= [X(t k ), IL M _(X k (t k ))] 



by the Poisson structure on M. that we will define in Section (3.3 1. This gives a hierarchy of 
commuting independent flows for k 



1.. 



1. From this we can sec that the trace of X is a 



Casimir with trivial flow. The solution of (3.3 1 can be found by factorization as in (2.25): Factor 

e t k x k (o) ^ 



(3.4) 



e ^(0) = n (t k )b(t k ) 



where n(t k ) is lower unipotent and b{t k ) is upper-triangular. Then 
(3.5) X(t k ) =n-\t k )X(0)n(t k ) 



solves (3.3) 



3.2.1. Characterization of blow-ups via Bruhat decomposition of G/B. Fix the eigenvalues Aj, and 
consider the level set A4\ consisting of all matrices in M. with spectrum (Ai, A„). In the case of 
distinct eigenvalues, Ercolani, Flaschka, and Haine |25j construct a minimal nonsingular compacti- 
fication of A4a on which the flows (3.3 1 extend to global holomorphic flows. The compactification 
is induced by an embedding of Ai\ into the flag manifold G/B, where G = SL(n,C) and B is the 
upper triangular subgroup of G. 



Proposition 3.1. (Kostant, [S3],) Consider the matrix 



( Ai 1 
A 2 1 



(3.6) 



£A 



\ 




••• ■■■ A„_i 1 
V A n ) 

in 6a. Every X 6 .Ma can be conjugated to e\ by a unique lower-triangular unipotent matrix L: 

(3.7) X = Le K L- 1 . 

This defines a map of A4\ into G/B: 

j A : Ma -» G/B 

X i ► L- 1 mod B . 



(3.8) 



This mapping is an embedding |53j . and the closure, Ja(-A4a); of its image is a nonsingular and 
minimal compactification of A^a |25j . Le t Lq be the unique lower unipotent matrix such that 
X(0) = Loe^L^ 1 . Then the solution (3.5 1 is conjugate to 6a as X(t k ) — n~ 1 (tk)L eALQ 1 n(tk), 



where L 1 n(t k ) is lower unipotent. Thus, X(t k ) is mapped into the flag manifold as 

(3.9) jA(X(t k )) = Lo l n{t k ) mod B 

(3.10) = L, 



mod B 



Notice that even at values of t k where (3.101 is not defined because the factorization (3.4 1 is not 



possible, the equivalent expression (3.10) is defined. In this way, the embedding of X(t k ) into G/B 



completes the flows through the blow-up times. This gives a compactification of A4a in G/B. [2"5] 
uses this embedding to study the nature of the blow-ups of X(t k ). 
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To illustrate this in a simple case, consider Example 1.1 from the Section [2T4] The isospectral 
set of 2 x 2 Hessenberg matrices with both eigenvalues zero is embedded into the flag manifold 
SL(2,C)/B, which has the cell decomposition 

(3.11) SL(2,C)/B= NB/B U N ^ B/B . 

Here N is the subgroup of lower unipotent matrices. The big cell, NB/B, contains the image of 
the flow X(t) whenever this flow is defined, that is, whenever the factorization e tx " — n(t)b(t) is 
possible. At t — — 1, where X(t) is undefined, the embedding j A completes the flow through the 
singularity. The image j\(X(t)) passes through the flag L^e^ 1 *^ at time t = — 1, which is the 



cell on the right in (3.11 1. 

The cell decomposition ( 3. 1 1[ > is a special case of the cell stratification of the flag manifold G/B 



known as the Bruhat decomposition. This decomposition is defined in terms of the Weyl group, W, 
as 

(3.12) G/B= (J NwB/B. 

In the present case of G = SL(n,C), W is essentially the group of permutation matrices. Thus, 
the Bruhat decomposition partitions flags according to which permutation matrix w is needed to 
perform the factorization X — nwb for X € G with n € N and b 6 B. At all values of tk for which 
the flow X(tfc) is defined, j\ sends Xfa) into the big cell of the Bruhat decomposition, since w is 



the identity. When the factorization (3.4 1 is not possible at time tk — t, e tx ^ can be factored as 



(3.13) e ix "^ = n(t)wb(t) 



for some permutation matrix w. In this case, the flow (3.101 enters the Bruhat cell NwB/B at time 
tk — i. [25] characterizes the Laurent expansion of each pole of X(t\) in terms of the Bruhat cell 
that the solution enters at the blow-up time. 



It is also seen in [25] that for k = 1, n — 1, the flows (3.10 1 generate a (C*) n torus action on 
G/B that has trivial isotropy group at every j\(X) with gk(X) ^ for all k. The orbit through any 
such point is open and dense in ja(AIa), and the closure of this orbit is the minimal compactification 
of j A (M A ) m G/B. 

3.2.2. Compactification of iso-level set with arbitrary spectrum. Here M and A^a are again defined as 
in Section |3~2] Shipman [75] uses a different embedding, referred to as the Jordan embedding, of AIa 
into G/B to describe the compactification of an isospectral set A4\ with arbitrary spectrum. The 
advantage of the Jordan embedding is that the maximal torus generated by the flows is diagonal if 
the eigenvalues are distinct and a product of a diagonal torus and a unipotent group when eigenvalues 
coincide. The orbits of these groups, specifically the torus component, are easily studied by taking 
their images under the moment map, as explained below. This leads to a simple description of 
the closure of A^a in terms of faces of the moment polytope. Recall that in the real tridiagonal 



Hessenberg form of the Toda lattice studied by Kodama and Ye [51] (see Section 2.4 1, the flows 
through an initial matrix X preserve the sign of each that is not zero and preserve the vanishing 
of each that is zero. The open subset of the isospectral set where no g\. vanishes is partitioned 
into 2 n_1 components, according to the signs of the gu- The compactification of the isospectral set 
is obtained by completing the flows through the blow-up times and pasting the 2 n_1 components 
together along the lower-dimensional pieces where one or more gu vanishes, producing a compact 
manifold. In contrast to this, when X is complex, A1a is no longer partitioned by signs of the g^\ 
there is only one maximal component where no g^ vanishes. The n — 1 flows through any initial X 
with g k ^ for all k generates the whole component, as was also observed in [2"5] . 
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To define the Jordan embedding, let Ca be the companion matrix of X, 



/0 








(3.14) C A = _ 

••• 1 

\s„ s„_i ••• s 2 0/ 

Here the s^-'s are the symmetric polynomials of the cigcvalucs Aj, i.e. 

n 

det(A7 - X) = A" - £ SjA"" J ' . 



Again, by [53], there exists a unique lower unipotent matrix L such that X = LCaL^ 1 . In particular, 
all elements of M\ are conjugate. Since the companion matrix has a single chain of generalized 
eigenvectors for each eigenvalue, any matrix in Jordan canonical form that is conjugate to it contains 
one block for each eigenvalue. 

Following |70j . fix an ordering of the eigenvalues, and let J be the corresponding Jordan matrix. 
Then C = \VJW~ 1 where W is a matrix whose columns are (generalized) eigenvectors of C, where 
each eigenvector has a 1 in the first nonzero entry and the vectors are ordered according to the 
chosen ordering of eigenvalues, with generalized eigenvectors ordered successively. Once W is fixed, 
for 1 e jUa, we can write X — LW JW~ X L~ X . The Jordan embedding is the mapping 

7A : Ma — » G/B 

(3.15) 

X i — > W^L- 1 modB 



That this is an embedding follows from the results in |53j 



Under this embedding, the flows X(t k ) = n (t k )X(0)n(t k ) in (3.51 with e k 



a t k x k -\Q) = 



as in (3.4 1 and X(0) — LWJW 1 L 1 generate a group action as follows: 



n(t k )b(t k ) 



lA{X{t k )) = W- 1 L- 1 n k {t k ) modB 

exp[ t k X k (0) } modB 



W- 1 L- 1 
W- 1 L- 1 



W- 1 L- 1 L Wexp[ t k Jl ] W~ 



exp[ t k (LWJ A W- 1 L- 1 ) k 



mod-B 



1 L- 1 modB 



The flows exp[ t k J^ ] for k = 1, ...,n - 
subgroup as A. j. 

Aj has r blocks along the diagonal 

{ ( s 



exp[ t h J% } W- 1 L- 1 modB 
n — 1 generate the centralizer of J in SL(n, < 



We denote this 



(3.16) 



A,, = { 



V 



s 



s e 



Xl, ■ 



J 



i = 1, . . . , r, where di is the multiplicity of the eigenvalue in the ith block of J. All the blocks together 
contain n — r independent entries in C above the diagonal and r entries in C* on the diagonal, where 
the product of the diagonal entries is 1. Aj is a semi-direct product of the diagonal torus Kj, 
obtained by setting all the entries above the diagonal equal to zero, and the unipotent group Uj, 
obtained by setting all the diagonal entries equal to 1. The subgroup of Aj that fixes every point 
in G/B is the (discrete) subgroup D of all constant multiples of the identity. The quotient Aj/D 
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has the manifold structure (but not the group structure) of (C*) r_1 x C' i_r . When r — n (distinct 
eigenvalues), Aj is the maximal diagonal torus. The compactification, Ma, of .Ma in G/B is the 
closure of one generic orbit of Aj. Its boundary is a union of non-maximal orbits of Aj. |72j uses 
the moment map of the maximal torus action, which sends SL(n,C)/B to a polytope in K" _1 , to 
identify each component of the boundary of Ma with a specified face of the polytope, as described 
above. 

First we describe the boundary of Ma in M. Let a be a subset of {<?i, <? n _i}, and denote 
by M\ the subset of Ma on which exactly the gi in a are zero. These subsets form a partition of 
Ma where the complex dimension of M% is equal to the number of gi that do not vanish. There is 
one maximal component, on which no gk vanish, and one component consisting of the fixed points, 
where all the vanish. 

Let X <G Mj[. The blocks on the diagonal of X where no gk vanish are full tridiagonal Hessenberg 
matrices of a smaller dimension (all entries on their first subdiagonals are nonzero). The union of the 
eigenvalues of these blocks, counted with multiplicity, is the spectrum A. Let P(A) be a partition of 
A into subsets A = Ai U ... U A q , where A& is the spectrum of the fcth block along the diagonal of X, 
and denote by M P t A \ the component of M% where A is partitioned among the blocks according to 



P(A). The Toda flows (3.5) through X preserve the spectrum of each block and therefore respect the 
partition Ma — U-Mp^. The moment map, described above, gives a one-to-one correspondence 
between the components Mp^ and particular faces of a certain polytope [72] , 

To see this, let Kj be the torus that lies along the diagonal of Aj. Kj is a subtorus of the 
maximal diagonal torus H . Its Lie algebra, JCj, is the kernel of a subset Aj of the simple roots; this 
determines the subset Sj C S of reflections in the hyperplanes perpendicular to the roots in Aj. Sj 
generates a subgroup Wj of W . The elements in W J = {w € W : l{sw) > l(w) Vs € Sj} are the 
coset representatives of minimum length in the quotient Wj\W. The following result is proved in 

m 

Proposition 3.2. [T2J The composition /io ^ : Ma ~* gives a one-to-one correspondence 
between the components M P ^ that partition Ma and the faces of the moment polytope with at 
least one vertex in W J . The complex dimension of the component is equal to the real dimension 
of the face. In particular, the maximal orbit in Ma corresponds to the full polytope, and the fixed 
points of Aj in Ma correspond to the vertices in W J . 



3.3. The full Kostant-Toda lattice. Here we consider full complex Hessenberg matrices 

/ * 1 ••• \ 



(3.17) 



X = 



1 



V 







/ 



The set of all such X is denoted e + where e is the matrix with l's on the superdiagonal and 
zeros elsewhere and i3_ is the set of lower triangular complex matrices. 

With respect to the symplectic structure on e+B- defined below, the Toda hierarchy (3.3) with X 



as in (3.171 turns out to be completely integrable on the generic leaves. The complete integrability 
is observed in [SB] by extending the results of [22] to e + B-. For n > 3, the eigenvalues of the 
initial matrix do not constitute enough integrals for complete integrability; a generic level set of 
the constants of motion is a subset of an isospectral set that is cut out by additional integrals and 



Casimirs, which can be computed by a chopping construction given in Proposition (3.3 1 
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Ikeda (40] studies the level sets in e + £>_ cut out by fixing only the eigenvalues. He finds a 
compactification of an isospectral set with distinct eigenvalues, showing that its cohomology ring is 
the same as that of the flag manifold SL(n,C)/B. This work differs from the previous works [HU] , 
[72] , and [5T] that compactify tridiagonal versions of the Toda lattice in that it does not use the 
Toda flows directly in producing the compactification. 

To describe the symplectic structure on e + B- , write 

(3.18) sl(n,C) =M- 

where A/1 and B + are the strictly lower triangular and the upper triangular subalgebras, respec- 
tively. With a non-degenerate inner product (A,B) = tr(AB) on sl(n, C), we have an isomorphism 
sl(n,£) = sl*(n,C) and 

sl*(n, C) = ATI © B% = Bi 8 ■ 

With the isomorphisms 

B+ = A/1: = B-, Afl = B+ = A/+ 

we identify 

e + z?_ = b; , 

which defines the phase space of the full Kostant-Toda lattice. On the space B* + , we define the 
Lie-Poisson structure (Kostant-Kirillov form); that is, for any functions f,h on B+, 

{f,h}(X) = (L,[U B+ Vf,U B+ Vh}), for XeB* + , 

where (Y,V/) = lim e _>o ^/(X+eY). This Lie-Poisson structure gives a stratification of B* + = e+B-. 
The stratification of the Poisson manifold B* + with this Lie-Poisson structure is complicated, having 
leaves of different types and different dimensions. 

Denote by B + the upper-triangular subgroup of SL(n, C) and by g ■ Y the adjoint action of 
SL(n,C) on sl(n, C). Then, through the identification of e + B- with B\, the abstract coadjoint 
action of B + on B\ becomes 

AdlX = e + Tle^- 1 -(X-e). 
The symplectic leaves in e + B- arc generated by the coadjoint orbits and additional Casimirs. 

In general, the dimension of a generic leaf is greater than 2(n — 1), and more integrals are needed 
for complete integrability. The chopping construction used in |22| to obtain a complete family of 
integrals for the full symmetric Toda lattice is adapted in [26] to find a complete family of integrals 
for the full asymmetric Toda lattice. 



Proposition 3.3. |26j Choose X G e + B-, and break it into blocks of the indicated sizes as 





k 


n-2k 


k 


k 1 




x 2 


x 3 


X = n-2k 


x 4 


x 5 


x e 


k 


\x 7 




x 9 



where k is an integer such that < k < [ ^ n 2 ]■ If d£&{X*?) ^ 0, define the matrix ipk(X) by 
(j>k(X) = X 5 - XiX^X s £ Gl(n - 2k, C), k ^ 0, 

MX) = x. 

The coefficients of the polynomial det(\—cf>k(X)) — A™ _2fe +/ifcA ,l_2fe_1 + - • -+I n -2k,k are constants of 
motion of the full Kostant-Toda lattice. The functions I\k are Casimirs on e + B-, and the functions 
Irk for r > 1 constitute a complete involutive family of integrals for the generic symplectic leaves of 
e + B- cut out by the Casimirs I\k. These integrals are known as the k-chop integrals. 
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The fc-chop integrals I r k are equivalent to the traces of the powers of <f>k{X). The Hamiltonian 
system generated by an integral I(X) is 

d 



(3.19) 



dt 



X = [X, ILv_(V/(X))] . 



When I(X) is one of the original Toda invariants Hk(X) = -j^jtr(X k + 1) (a 0-chop integral), the 
flow is 



jX = [X, U u _{X k )} . 



(3.20) 

The solution may again be found via factorization |26] . Let 

e tVI(X ) = n{t)b{t) > 

with n(t) and b(t) lower unipotent and upper-triangular, respectively. Then 

X(t) = rC l {t)Xan{t) . 



Let X belong to (e + S_)a- Recall from Section |3.2| that there exists a unique lower unipotent 
itrix 

(3.21) 



matrix L such that X = LC\L , where C is the companion matrix (3.14). The mapping 

Ma — ► SL(n,C)/B 
X i — ► L- 1 modB 



ca 



is an embedding |53j . referred to as the companion embedding. Its image is open and dense in the 
flag manifold. Under this embedding, the n — 1 flows of the 0-chop integrals j:trX k generate the 
action of the centralizer of Ca in SL(n,C) (the group acts by multiplication on the left). 
When the Xi are distinct, Ca = VAV -1 , where V is a Vandermonde matrix, and 



X = LVAV~ X L~ 



The embedding 
(3.22) 



tf A : (e + B_) A 
X 



SL(n,C)/B 
V^L- 1 modB 



is a specific cas e of t he Jordan embedding (|3.15|) when the eigenvalues are distinct. In this case, the 

^tx{X k ) for k = l,.. 



1 is the 



group Aj (see (3.16 1) generated by Hamiltonian flows of Hk 
maximal diagonal torus. <J a is therefore referred to as the torus embedding. 

When the values of the integrals are sufficiently generic (in particular, when the eigenvalues of 
each fc-chop are distinct), Ercolani, Flaschka, and Singer show how the flows of the fc-chop 
integrals can be organized in the flag manifold by the torus embedding. (The companion embedding 
gives a similar structure, but the torus embedding is more convenient since the group action is 
diagonal.) The guiding idea in [26] is that the fc-chop integrals for SL(n,C) are equivalent to the 
1-chop integrals for SL(n — 2{k — 1),C). Let SL(m,C)/B denote the quotient of SL(m,C) by its 
upper triangular subgroup, and let SL(m,C)/P denote the quotient of SL(m,C) by the parabolic 
subgroup P of SL(m, C) whose entries below the diagonal in the first column and to the left of the 
diagonal in the last row are zero: 



P 



( * 






* 


\ 









* 











* 




V o 


• 


• 


* 


/ 
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The 1-chop integrals I r \ depend only on the partial flag manifold SL(n,C) / P. In this partial flag 
manifold, a level set of the 1-chop integrals is generated by the flows of the 0-chop torus. The 1-chop 
flows generate a torus action along the fiber of the projection 

SL(n-2,C)/B -» SL(n,C)/B 

I 

SL(n,C)/P . 

In this fiber, the 2-chop integrals depend only on the partial flag manifold SL(n — 2, C) /P, where 
a level set of the 2-chop integrals is generated by the 1-chop torus. This picture extends to all the 
fc-chop flows. [26 builds a tower of fibrations 

SL(n-2(k+l),C)/B -> SL(n-2k,C)/B 

i 

SL(n-2k,C)/P 

where the fc-chop flows generate a level set of the (k + l)-chop integrals in the partial flag manifold 
SL(n — 2k, C)/P and the (k + l)-flows act as a torus action along the fiber, SL(n — 2(k + l),C)/B. 
In the end, the closure of a level set of all the fc-chop integrals in SL(n, C)/B is realized as a product 
of closures of generic torus orbits in the product of partial flag manifolds 

(3.23) SL(n, C)/P X SL(n - 2, C)/P X • • • X SL(n - 2M, C) jP 

where M is largest k for which there are fc-chop integrals. 

In |33] , Gekhtman and Shapiro generalize the full Kostant-Toda flows and the fc-chop construction 
of the integrals in Proposition |3.3| to arbitrary simple Lie algebras, showing that the Toda flows on 
a generic coadjoint orbit in a simple Lie algebra g are completely integrable. A key observation in 
making this extension is that the 1-chop matrix 4>i(X) can be obtained as the middle (n — 2) x (n — 2) 
block of Ad-p(x)(X), where T(X) is a special element of the Borel subgroup of G. This allows the 
authors to use the adjoint action of a Borel subgroup, followed by a projection onto a subalgebra, 
to define the appropriate analog of the 1-chop matrix. 

Finally, we note that full Kostant-Toda lattice has a symmetry of order two induced by the 
nontrivial automorphism of the Dynkin diagram of the Lie algebra sl(n, C). In terms of the matrices 
in e + B-, the involution is reflection along the anti-diagonal. It is shown in [67] that this involution 
preserves all the fc-chop integrals and thus defines an involution on each level set of the constants 
of motion. In the flag manifold, the symmetry interchanges the two fixed points of the torus action 



that correspond to antipodal vertices of the moment poly tope under the moment map (3.1 1 



Example 3.1. In this example, we demonstrate the complexity of the Poisson stratification of e+B- 
for n = 3 and n = 4. The table of symplectic leaves of all dimensions has been calculated in notes by 
Stephanie Singer, a co-author of |26j . as given below. On the leaves of lower dimensions, the fc-chop 
integrals are dependent. 
When n = 3, 

( / t i n \ 

B - = { i m h 1 j : E/i = 

i=l 

Its symplectic leaves are listed in the following table. The Casimirs are constants of motion that 
generate trivial Hamiltonian flows. The value of each Casimir is fixed on a given symplectic leaf. 
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h = 
h, gi = 
h, ,g 2 = 
h, 5i, 52 = 



5i52 7^ 
52 ^0 
51^0 



9192 
h 



fl 
h 

h, h, C/3) 



4 
4 
2 
2 




E/, = o 

i=l 



On each 4-dimensional leaf, the functions ^tr(X 2 ) and |tr(X 3 ) provide a complete family of integrals 
(Hamiltonians) for the Toda hierarchy. 

When n = 4, the symplectic stratification is already much more complicated. Here, 

{ ( h i o o \ 

5i /2 1 
hi 52 h 1 

k V k h 2 93 U J 

The table of symplectic leaves is as follows: 



9'jhi 
k 



Dimension 



Closed Conditions 




Open Conditions 


Casimirs 






fc(fc 52 - h x h 2 ) 56 


hi = -h - h + Sih ^ 


fc.92 — hifi2 = 




fc^O 


Ci = /2- 2 ^,C 2 = / 3 - 


fc = 




hih 2 {g\h 2 + 53/11) 7^ 




fc, /ii = 






r — f 9293 

C 3 - 73 - -S^- 


fc, /l 2 = 




/ii53 + 


C4 = /2 - 


fc, 5i /i 2 + 53^1 = 




/ll^2 7^ 


C 5 = /1 + /s - 9 ~ft 


fc, /ii, h 2 = 




515253 7^ 




fc, /ii, 51 = 






/i,C 6 = / 3 -2|f 


fc, fti, 51, h 2 = 




5253 7^ 


/1 


fc, h 2 , 53 = 




hi 7^0 


.A, c 7 = / 2 - 


fc, /l2, 53, /ii = 




5i52 7^ 


/4 


fc, /ii, h 2 , 52 = 




5i53 7^ 


/l + /2 (= -(/3 + .A)) 


fc, /li, /i 2 , 51, 52 = 


= 


53 7^0 


fl, h 


fc, /ii, h 2 , 52, 53 = 


= 


517^0 


f3, fi 


fc, hi, h 2 , 51, 53 = 


= 


52 7^0 


fl, fi 


fc, fti, ^2, 5i, 52, 


53 = 




fl, j*2, /3, (A) 



On the maximal leaves, of dimension 8, the functions | trX* for fc = 2, 3, 4 provide three constants 
of motion. One 1-chop integral is needed to complete the family. 
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1120 




1102 

Figure 3.2. The moment polytope A4 for SL(4,C)/P. The vertices represent the 
weights Li~Lj which are expressed by iiLi-H2-£2 + *3£3+*4£4 using Li+- • - + £4 = 
0, e.g. 2110 means L\ — L4. 



3.4. Nongeneric flows in the full Kostant-Toda lattice. When eigenvalues of the initial matrix 
in e + B- coincide, the torus embedding (3.22) is not defined since any matrix in e + B- has one 
Jordan block for each eigenvalue. In the most degenerate case of non-distinct eigenvalues, that is, 
when all eigenvalues are zero, the isospectral set can be embedded into the flag manifold by the 
companion embedding (3.211. Under this embedding, the 0-chop integrals generate the action of 
the exponential of an abelian nilpotent algebra |65j . The 1-chop integrals are again defined only in 
terms of the partial flag manifold SL(n, C) /P. Fixing the values of each 1-chop integral produces a 
variety in the flag manifold. The common intersection of all these varieties turns out to be invariant 
under the action of the diagonal torus and has a simple description in terms of the moment polytope 

[70] considers level sets where the eigenvalues of each 4>k(X) are distinct but one or more eigen- 
values of 4>j{X) and <fij+i(X) coincide for one or more values of j. In this situation, the torus 



orbits generated by the fc-chop integrals in the product (3.23 idegenerate into unions of nongeneric 
orbits. The nature of this splitting can be seen in terms of the moment polytopes of the partial flag 



manifolds in (3.231 



Recall the definition of the moment map // in (3.1). Here G is SL(n,C) and V is the adjoint 
representation. V may be realized as the subspace of C" ® A"~ 1 C" with e i ® e* = 0, where {e^} 
is the standard basis of C, and e* = (— l) J+1 ei A ... A e,_i A e i+i A . . . A e„. The partial flag manifold 
SL(n, C) jP is the orbit of G through [ei ® e* ] in ¥(V). The weight of ® e* is Li — Lj, where L% is 
the linear function in Ti.* that sends an element of TL to its fcth diagonal entry. The weights Li — Lj 
with i j are the vertices of the weight polytope of V, which we denote by A„. These vertices are 
the images under the moment map of the fixed points of the complex diagonal torus. The image of 
the closure of a torus orbit under moment map is the convex hull of the weights corresponding to 
the fixed points of the torus in the closure of the orbit. The real dimension of the image is equal to 



the complex dimension of the orbit [5]- Figure 3.2 shows the example of the moment polytope A4. 

An element gB in SL(n,C)/P represents the partial flag V 1 C V™ -1 C C™ where V 1 is the span 
of the first column of g and U n_1 is the span of the first n — 1 columns. There are two natural 
projections from SL{n, C)/P to the projective space CP™ -1 and its dual (CP" -1 )* that send a 
partial flag to the line V 1 and to the hyperplane V™ -1 , respectively. Let iTi and tt* be projective 
coordinates on CP" -1 and (CP" -1 )*. 

The coordinates "k^-l- that come from the embedding of SL(n,C) / P into ¥(V) by the moment 
map (3.1 1 are projectively equal to the products 7Tj7r* for i ^ j: [^Li-L^ijtj = [^i^jhjij ■ 
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Figure 3.3. The polytopes A 4 (\fc) and A 4 (\k*), obtained by splitting A 4 along 
an interior hexagon. Each is missing the vertices of one triangular face. 

At each fixed point of the diagonal torus in SL(n,C)/P, exactly one -Ki and one tt* does not 
vanish. Those where ttu 7^ correspond to the vertices — Li with i ^ k, whose convex hull is an 
(n — 2)-dimcnsional face of A„, which we denote as A n (k). The fixed points where n%. 7^ correspond 
to the vertices L, ; — of the antipodal face, A„(fc*). The polytope of an (n — l)-dimensional torus 
orbit where 7Tfe or 7r£ is the only vanishing coordinate is the convex hull of the vertices remaining 
after the vertices of the face A„(fc), respectively A n (fe*) are removed. These polytopes are denoted 
A„(\fc) and A n (\k*), respectively. They are congruent polytopes, obtained by splitting A„ along 
the hyperplane through the vertices Li — Lj with i, j =/= k. The convex hull of these vertices is an 
(n — 2)-dimcnsional polytope in the interior of A„, which we denote as A„(\fc\fc*). We will refer 



to the pair A n (\fc) and A n (\fc*) as a split polytope. In Figure 3.3 we illustrate the example of the 
split polytope A 4 (\fc) and A 4 (\fc*) [70]. When two or more such splittings occur simultaneously, 
the collection of resulting polytopes will also be called a split polytope. 

Proposition 3.4. |70j Let T be a variety in SL(n,C)/P defined by fixing the values of the 1-chop 
integrals I r \, including the Casimir, where the values are chosen so that exactly one eigenvalue, say 
Aio, of X is also an eigenvalue of 4>\(X). Then T is the union of the closures of two nongeneric 
torus orbits, O l and O l , on which -Ki, respectively n*, is the only coordinate that vanishes. The 
images of their closures under the moment map, A n (\z) and A n (\i*), are obtained by splitting A n 
along the interior (n — 2) -dimensional face A n (\i\i*). When exactly p eigenvalues of X are also 
eigenvalues of (j)\{X) (p < n — 2), then T is the union of the closures of 2 P nongeneric (n — 1)- 
dimensional orbits whose images under the moment map are the polytopes obtained by splitting A n 
simultaneously along p interior faces A n (\j\j*). 

This result extends to the fc-chop flows as follows: 

Proposition 3.5. jTO] If p eigenvalues of (jj^^X) and (f>k-i(X) coincide, then the generic orbit of the 



diagonal torus that generates the {k— 1) -chop flows in the component SL(n — 2(k—l), < C)/P of (3.23) 
becomes a union of2 p nongeneric orbits. Since the moment map on the product \ 3.23\ is the product 



of the component moment maps, the moment map on the product of partial flag manifolds takes a 



level set in (3.23) to a product of full and/or split polytopes, depending on where the coincidences of 



eigenvalues occur. 

When a level set of the constants of motion is split into two or more nongeneric torus orbits, there 
are separatrices in the Toda flows that generate the torus action. The faces along which the polytope 
is split are the images under the moment map of lower-dimensional torus orbits (the separatrices) 
that form the interface between the nongeneric orbits of maximum dimension. The flow through an 
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initial condition in one maximal orbit is confined to that orbit. It is separated from the flows in the 
complementary nongeneric orbits by the separatrices. |68j determines the monodromy around these 
singular level sets in the fiber bundle of level sets where the spectrum of the initial matrix is fixed 
with distinct eigenvalues and the remaining constant of motion / (equivalent to the determinant of 
the 1-chop matrix) is allowed to vary. The flow generated by / produces a C*-bundle with singular 
fibers over the values of /. The singularities occur both at values of / where an eigenvalue of the 
1-chop matrix coincides with an eigenvalue of the original matrix and at values of / where the two 
eigenvalues of the 1-chop matrix coincide. In a neighborhood of a singular fiber of the first kind, the 
monodromy is characterized by a single twist of the noncompact cycle around the cylinder C* . Near 
a singular fiber of the second kind, the monodromy produces two twists of the noncompact cycle. 
This double twist is seen in the simplest case when n = 2 near the fiber where the two eigenvalues 
of the original matrix coincide; it as described in detail in [66 . 



When eigenvalues of <j)\{X) coincide, the torus embedding (3.221 generalizes to the Jordan embed- 



ding (3.151, under which the 1-chop flows generate the action of the group Aj in (2.201, a product 
of a diagonal torus and a nilpotent group. The general structure of a level set of the 1-chop integrals 
with this type of singularity is not known, in part because the orbit structure of Aj in the flag 
manifold is not understood in sufficient detail. When the eigenvalues are distinct, Aj is a diagonal 
torus, and the closures of its orbits are toric varieties [SS]- The structure of torus orbits in flag 
manifolds is well- understood; see for example [5], [5U], and [53]. The closures of orbits of Aj in the 
flag manifold are generalizations of toric varieties, and much less is known about them. 

The fixed points of the actions of the groups Aj are studied in [69], and the fixed point sets of the 
torus on the diagonal of Aj are characterized in [71]. If Aj has r blocks along the diagonal, where 
the dimension of the ith block is dj, then the maximal diagonal subgroup of A has gcd(ii 1 , . . . ,d r ) 
connected components [75 . The subgroup of Aj that fixes all points in the flag manifold is the 
discrete group D consisting of constant multiples of the identity where the constants are the nth 
roots of unity; the group Aj/D then acts effectively on the flag manifold. [73] describes the fixed- 
point set of the unipotent part of Aj, giving an explicit way to express it in terms of canonical 
coordinates in each Bruhat cell. In the case where all eigenvalues coincide, Aj is equal to its 
unipotent part. [73] shows that the action of the group in this case preserves each Bruhat cell and 
that its orbits in a given cell are characterized by the "gap sequence" of the permutation associated 
to the cell. 

4. Other Extensions of the Toda Lattice 

4.1. Isospectral deformation of a general matrix. In the full Hessenberg form of the Toda 
lattice, the matrix is diagonalizable if and only if the eigenvalues are distinct. Kodama and Ye 
generalize this in |52j , where they consider an iso-spectral deformation of an arbitrary diagonalizable 
matrix L. The evolution equation is 

(4.1) j t L=[P,L]; 
P is defined by 

(4.2) P = U(L) = (L) >0 -(L) <0 , 

where (£)>o(<o) is the strictly upper (lower) triangular part of L. [52] establishes the complete 



integrability of (4.1 1 using inverse scattering, generalizing the method used in [37] to solve the full 



symmetric real Toda lattice. The method yields an explicit solution to the initial-value problem. 



The general context of the flow (4.1 1 includes as special cases the Toda lattices on other classical Lie 
algebras in addition to sl(n, K), which is most closely associated with Toda's original system. In this 
regard, Bogoyavlensky in [J2] formulated the Toda lattice on the real split semisimple Lie algebras, 
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which are defined as follows (the formulation below is in the Hessenberg (or Kostant) form): Let 
{h a . , e± ai : i = 1, ...,/} be the Chevalley basis of the algebra g of rank I, i.e. 

where (CV/)i<ij<z is the Cartan matrix and CV, = ct;(/i Qj . ). Then the (non-periodic) Toda lattice 
associated with the Lie algebra g is defined by the Lax equation 



(4.3) 



dL 
~dt 



[A,L], 



where L is a Jacobi element of g and A is the A/"_-projection of L, 



i(*) = E fi(t)h ai + E(9i(t)e 



A(t) = —Hjsf_ L(t) = - £ 9i{t) e_ Ql . 

i=l 

The complete intcgrability is based on the existence of the Chevalley invariants of the algebra, and 
the geometry of the isospectral variety has been discussed in terms of the representation theory of 
Lie groups by Kostant in |54j for the cases where gi are real positive, or complex. The general case 
for real gi's is studied by Casian and Kodama jT21[TB], which extends the results in the sl(n, K.) Toda 
lattice in the Hessenberg form (see Section 2.4 1 to the Toda lattice for any real split semisimple Lie 
algebra. 



The Lax equation (4.3 1 then gives 



dt 

dgi 
dt 




from which the r-functions arc defined as 



(4.4) 



f k (t) = ~]n Tk (t), 



i 



-Gki 



In the case of g = sl(n,M), those equations are (2.281 and (2.291 (note here that the superdiagonal of 
L(t) is diag(/i — /2, — /3, •••,// — /z+i) with n = I + 1). Those extensions have been discussed by 
many authors (see for example [36 , 60 ) . One should note that Bogoyavlensky in [T^] also formulates 
those Toda lattices for affine Kac-Moody Lie algebras, and they give the periodic Toda lattice. There 
has been much important progress on the periodic Toda lattices, but we will not cover the subject 
in this paper (see for example [21 |3j [24j [621 163]). 

From the viewpoint of Lie theory, the underlying structure of the integrable systems is based on 
the Lie algebra splitting, e.g. sl(n) — £?_ © so(n) (the QR-decomposition) for the symmetric Toda 
lattice, and sl(n) = B+ ®J\f- (the Gauss decomposition) for the Hessenberg form of Toda lattice. 
Then one can also consider the following form of the evolution equation, 



(4.5) 



d 
dt 



L=[Q k ,L] with Q k = U 3l (L k ). 



where 2i is a subalgebra in the Lie algebra splitting sl(n) — Qi ©02- In this regard, we mention here 
the following two interesting systems directly connecting to the Toda lattice: 
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(a) The Kac-van Moerbeke system [42]: We take g\ = so(2n), and consider the equation for L 6 
so(2n) (recall that L is a symmetric matrix for the symmetric Toda lattice) . Since L 2fe_1 G so(2n), 
the even flows are all trivial. Let L be given by a tridiagonal form, 



L = 



( o 

-ax 


V o 



Oil 







Oil 



\ 





"C*2n-1 



Then the even flows are the Kac-van Moerbeke hierarchy. 
Skew^ 2 - 7 )), where the first member of i2-flow gives 

9a k , 2 2 



at 2 



- a k+i) > 



tt2n-l 

/ 



i. 



G so(2n,R) 



[IL (L 2i ), L] (recall that n so (L 2 ^) 



,2n-l 



with ao = a 2rl = 0. This system is equivalent to the symmetric Toda lattice which can be written 
as (4.5) for the square 1? . Note here that 1? is a symmetric matrix given by 



r (i) 



1 




T (2) 




1 



where rW, for i = 1,2, are n x n symmetric tridiagonal matrices given by 



rp(i) 



(bf 

.« 

1 





V o 





(i) 



o \ 




w 

n-1 

(i) 
'n-1 



1.0 



with a 



(i) 



*2fe-2 



(2) 



l 2fc-l> "fc 



l(2) 



*2fc-l 



«2fe (see [35]). 



aafcOafc+ij and 

Then one can show that each T^> gives the symmetric Toda lattice, that is, the Kac-van Moerbeke 
hierarchy for L 2 matrix splits into two Toda lattices, 



m 



2j 



= [n so (T«)^',T«] i= 1,2. 



The equations for TW are connected by the Miura-type transformation, with the functions [al , &L ), 
through the Kac-van Moerbeke variables (see [35] )• 

(b) The Pfaff lattice for a symplectic matrix [TJ [48] : The Pfaff lattice is defined in the same form 
with gi = sp(2n) and L in the Hessenberg form with 2x2 block structure. In particular, we consider 
the case L € sp(2n) having the form, 



L = 



/0 si 
h 



ax 




02 \ 



ax 


s 2 
b 2 




o 2 










V ° 2 


o 2 




s„ 
6 n / 



£ sp(2r 
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where O2 is the 2x2 zero matrix. The variables (dk,bk) and = ±1 are those in the indefinite 
Toda lattice. It should be noted again that the odd members are trivial (since L 2k ~ x £ sp(2n)), and 
the even members give the indefinite Toda lattice hierarchy |49j . Here one should note that L 2 can 
be written as 

,0 o) +i ® o 1, 



L 2 = L T 



where L is given by (2.331. Then one can show that the generator Q 2 j of the Lax equation is given 
by 



B 







B, 



1 



where B } = ^[(V) >0 - [V) <0 ]. Then the hierarchy f f L = [Q 2j ,L) gives the indefinite Toda lattice 
hierarchy. 

4.2. Gradient formulation of Toda flows. In [7 , Bloch observed that the symmetric tridiagonal 



Toda equations (2.6) can also be written in the double-bracket form 

d 



(4.6) 



dt 



L(t) = [L(t), [L(t), N]] 



where N is the constant matrix diag(l, 2, n) and L is as in (2.7l. He showed that this double- 
bracket equation is the gradient flow of the function f(L) = tr(LN) with respect to the normal 
metric on an adjoint orbit of SO(n). The normal metric is defined as follows: Let k( ,) = —(,) be 
the Killing form of a scmisimplc Lie algebra g, and decompose g orthogonally relative to ( , ) into 
g = g L gL where g^ is the centralizcr of L and g L = Im &d(L). For X £ g, denote by X L the 
projection of X onto g L . Then given two tangent vectors to the orbit at L, [L,X] and \ L, Y] , the 
normal metric is defined by ([L, X], [L,Y])jf — (X L ,Y L ). Then the right hand side of (4.6 1 can 
be written as gradiJ = [L, [L,N]] for the Hamiltonian function H(L) — k(L,N) (Proposition 1.4 
in [S]). Thus the Toda lattice (2.6 1 is both Hamiltonian and a gradient flow on the isospectral set. 



Brockett shows in 1 13] that any symmetric matrix L(0) can be diagonalized by the flow ( |4.6| , and the 
flow can be used to solve various combinatorial optimization problems such an linear programming 
problems (see [37] for the connections of the Toda lattice with several optimization problems) . 
The flow (4.6 1 is extended in [5] and [S] to show that the generalized tridiagonal symmetric 



Toda lattice can also be expressed as a gradient flow. In Section 4.1 we give the equations of the 
generalized tridiagonal Toda lattice in the Hcsscnbcrg form on a real split semisimplc Lie algebra. 
The symmetric version of this is as follows (see [9]). Let g be a complex semisimple Lie algebra of 
rank I with normal real form g n . Choose a Chevalley basis {h att , e± ak : k = 1, . . . , 1} as in Section 



|4.1| The generalized tridiagonal symmetric Toda lattice is defined by the Lax equation 

d 



dt 



L(t) = [A(t),L(t)} 



where 



L(t)= £ b k (t)h ak 

k=l 



E ak(t)(e ak 
fc=i 



M~t) = E a k (t)(e ak -e_„J. 
fc=i 

This flow defines a completely intcgrable Hamiltonian system on the coadjoint orbit of the lower 

1 

Borel subalgebra of g n through E ( e a k ~ e -a k )- The Hamiltonian is H(L) — \K{L, L), where K is 

fe=i 

the Killing form. 
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[9] shows that (4.2 1 is a gradient flow with respect to the normal metric on the orbit. The gradient 
formulation in is given in the context of the compact form g u of g (see also the survey in 
Their key result is the following. 

Proposition 4.1. The gradient vector field of the function f(L) — K(LN) on the adjoint orbit in 
g u containing the initial condition Lq, with respect to the normal metric, is 



(4.7) 



~L(t) = [L(t), [L(t), N}} 



Now let H u be a maximal abelian subalgebra of g u , and take H = H u © iH u as the Cartan 
subalgebra of g. Choose a Chevalley basis for g as above. Bloch, Brockett, and Ratiu jS] show the 
following. 

Theorem 4.1. Let N be i times the sum of the simple coweights of g, and let 

i i 
L(t) = J2 ibk(t)h ak + J2 ia k (t)(e ak + e_ Q J . 
fc=i fe=i 



Then the gradient vector field (-i-D gives the flow of the generalized tridiagonal symmetric Toda 



lattice on the adjoint orbit in g u containing the initial condition Lq. Explicitly, 



N 



J2 ixkK 

k=l 



where (x\,...,x{) is the unique solution of the system 



Y, x k a p (h ak ) = -l,p= 
fe=l 

A list of the coefficients (xi, ...,£;) for all the semi-simple Lie algebras is given on p. 62 of [5]- 

to the generalized full symmetric Toda 



Proposition (14:. lh and Theorem (4.1l are extended in 



lattice and in |101 ITT] to the generalized signed Toda lattice to show that these extensions of the 
Toda flows are also gradient flows. 



5. Connections with the KP equation 

Here we give a brief review of the paper jS] whose main result is to show that the r-functions of 
the Toda hierarchy (2.44 ) with a symmetric tridiagonal matrix provide a new class of solutions of the 



Kadomtsev-Petviashvili (KP) equation. We also provide a geometric description of the r-functions 
in terms of the Grassmann manifolds (see |46j ) . 



5.1. The r-functions for the symmetric Toda lattice hierarchy. We return to the Toda lattice 



hierarchy (2.441 with symmetric tridiagonal matrix. The solution L(t) can be explicitly expressed 
in terms of the r-functions (2.131: Let us summarize the process of solution method based on the 



Gram-Schmidt orthogonalization. First we consider 

. 9 (t)=cxpQ#(L(0),t) 



where 6»(A,t) := £ t k X k with the fc-th flow parameter t k of the Toda hierarchy, i.e. 
fe=i 

BT 1 

— = [B k ,L] with B k = -Skcw(i fc ) . 

Ot k A 
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(Note here that we rescale the time tk — > tk/2.) Then define the matrix 
M(t) := g T (t)g(t) = 

= $(Q)e e ( A ' t )$ r (O) = f(^0?e 9 ( A ' t )) N ) 

V / l<i,j<n 

where $(0) = (<£°( A j))i<i,j<n is the eigenmatrix of L(Q), i.e. L(0)$(0) = $(0)A, and $(0) G 50(»). 
Since L(0) is a tridiagonal matrix, the entries mjj(t) := (<^</>°e e ( A ' t - ) ) can be written in terms of the 
moment by the Gram-Schmidt orthogonalization process (see [47] for the details), 

mi ,(t) = (\ i+j ~ 2 p(\) e°^) = t K +J ^ Pk e Sfe(t) , 

fc=i 

where p(X) = </>i(A) 2 with pk = p(Xk), and 0fc(t) = 9(Xk,t). In particular, we have 



(5.1) 



ri (t)=(p(A)e^))=E ft^'M 
fc=i 



Then the r-functions are given by the Wronskian of the set of functions of ri(t) and its x-derivatives, 



r fe (t)=Wr(Ti(t),r((t) 



.( fc -!)^ 



for k = 1, 2, . . . , n — 1 . 



Using the Binet-Cauchy theorem, one can write in the form 
(5.2) 



r /c(t) = E - - - , *fe)(t) , 

l<2i <---<2fc <n 



where 



(5.3) 



...,i fc ) = 



n (^-^) 
i<i<j<*-- 



[ £?(*!,...,<*) =Wr(e^,e°'» > ...,e e '*)- 
Here we assume the ordering in Aj as 

(5.4) Ai < A 2 < • • • < A„ , 

so that all £(«i, . . . , «&) are positive. As we show below that this form of has interesting geometric 
interpretation, which makes a connection with the KP equation. 



Example 5.1. In Figure 5.1 we illustrate the solutions (eti(t), et2(t), a^{t)) for the sZ(4, R) Toda 
lattice. The r-functions are given by 

Tl (t) = e ^ + • ■ ■ + eW\ 

r 2 (t)=E(A J '-A i )V*W+^( t ), 

T3 (t)= Y, [(A t -A J )(A J -A fe )(A l -A fc )]V-( t )+^( t )+ 9 '=( t ), 
with 0fc(t) = Afcii + A^2 + ^fc^3- The solutions afe(t) are then given by 



Tfc-l(t)T fc+ i(t) 

4{t) ■ 



4(t) 

which can be also written as a 2 = d 2 In Tk /dt\ . 

Each line in Figure 5.1 is given by the equation #fc(t) = 0,-(t) for some k and j, for example, in 
the middle graphs, two lines are 
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Figure 5.1. The solutions (ai, a 2 , 03) for sZ(4, K) Toda equation. The graphs show 
the contour lines of the solutions Ofc(ti, ta, ^3) f° r ti-t2 plane with £3 = 0. The left 
graph is for a%, the middle one for a 2 and the right one for 03. Each line is given 
by the balance 9j(t) — 9k(t) where 6j(t) = \jt\ + X 2 t 2 + A|t3. Here the eigenvalues 
(Ai, . . . , A4) are given by (—3, 0, 1, 2). For example, the lines in the middle graph 
are given by t\ + 2t 2 = and t\ — 2t 2 = 0. Those graphs also show soliton solutions 
of the KP equation (see Section |5.2| . 



(i) 0i(t) = 9 3 (t) which gives 

h + (Ax + A 3 ) t a + (A? + AiAg + \\) t 3 = , 

(ii) 2 (t) = 04 (t) which gives 

h + (A 2 + A 4 ) t 2 + (Aa + A 2 A 4 + \\) t 3 = . 

Each line indicates a balance between two exponential terms in the r- function, for example, 9j = 9^ 
shows the balance between e 6j and e 6k in n, or e 9j+ 1 and e k+0i in r 2 for some I. Then near each 
line for 9j — 9k, we have, from n ~ e ej + e 9fc , 

«*(t) = J(A fc - A.fsech 2 i(^(t) - <? fc (t)) . 

This is a soliton solution of the KP equation, and we will discuss some details of the connection to 
the KP equation in the next section. In this regard, the graphs in Figure [5~T| shows several solutions 
of the KP equation, and in particular those solutions indicate the soliton resonances (see [BJ USj |5T] 
for more details). 

5.2. The KP equation and the r-function. The KP equation is given by the following partial 
differential equation, 

(5 5) ® ( 4^" + ^ U + 12u^ uS ] + 3^ 2 " — 

dx \ dt dx 3 dx J dy 2 

where u = u(x, y, t) with {x, y) represents a coordinate of two space dimensions and t the time. The 
KP equation is one of the prototypical (2+l)-dimensional integrable equations, originally derived in 
43J as a model for small-amplitude, quasi two-dimensional waves in a weakly dispersive medium. It 
arises in many different physical applications including shallow water waves and ion acoustic waves 
in plasmas (for a review, see e.g. [4Tj). 

The goal of this section is to show that a large class of the solutions of the KP equation can be 
obtained in the Wronskian forms which contain the r-functions of the symmetric Toda lattice (this 
is based on [64], and see also [57] and Appendix in [6]). We start with a brief explanation of the 
inverse scattering scheme for the KP equation: Let C be a pseudo-differential operator defined by 

C = d + u 2 d- 1 +u 3 d- 2 + ••• , 
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where d is a derivation satisfying dd 1 = d 1 d = 1 and the generalized Leibnitz rule, 

oo / \ f)k f 

^ = £,(*) a? fOT 

(Note that the series terminates if v is a positive integer.) Then the following infinite set of equations 
is called the KP hierarchy: 

dC 

(5.6) — = [B n ,C], with B,- :=(/?)+, j = l,2,.... 

Here (£")+ represents the projection of £ n onto the polynomial (differential) part in d. For examples, 
the first three members of Bj are given by 

B 1 = 9, B 2 = d 2 + 2u 2 

B 3 = d 3 + 3u 2 d + 3(u 2 , x + u 3 ) 

where u 2tX — du 2 /dx. The equation B\ = d implies 

dC dC 

from which we identify t\ = x. The compatibility among the equations in (5.6 1 is given by the 
Zakharov-Shabat (ZS) equations, 

which are a direct consequence of the definition of Bj. In particular, the ZS equation with i = 2 and 
j = 3 gives 

2 <9w3 _ dua d 2 u 2 _ Q 
9x dt 2 dx 2 



( du 3 d 2 u 3 \ du 2 , _ <9 2 u 2 <9 3 u 2 , _ du 2 
— 2 V- .1 V- riu.n 



V dt 2 dx 2 J dt 3 ~^~^dt 2 dx dx 3 ~^^ U2 g x 



Eliminating u 3 from those equations, we obtain the KP equation (5.5) with the identifications u = u 2 
and y — t 2 ,t = t 3 . This means that any solution of the KP hierarchy is also a solution of the KP 
equation. 

Now writing C in the dressing form, 

(5.7) C = WdW~ 1 , with W = l+w 1 d- 1 +w 2 d- 2 + ■■■ , 
the KP hierarchy becomes 

dW 

(5.8) ■ w -=B j W-Wd } . 

atj 



Using (5.7l, the variables Ui can be expressed in terms of the Wj's; for example 

dw\ dw 2 dwi 

OX 

and so on. The equations for utj are, for example 



u 2 = --^— , u 3 = -— hu>i- , 

ox ox ox 



dwi _ ^ w , 9 2 wi ^ dw 2 
dt 2 1 dx dx 2 dx 

dw 2 _ ^ w dw i + ^ 2w 2 + 2 <9w 3 
dt 2 " dx dx 2 dx 



and so on. 
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To find some exact solutions, we consider a finite truncation of W, for some k > 1, 

W fe := 1 + w^- 1 + ■ ■ ■ + w k d- k . 



For example, the W-equation (5.8 I for the k — 1 truncation, i.e. W\ = l + W\d 1 , is just the Burgers 
equation, 



(5.9) 



-2W! 



dwi dwi 
dx dx 2 



which can be solved by a Cole-Hopf transformation, w\ = —d\nf/dx, leading to a linear diffusion 
equation, f t% = f xx . Note here that the Cole-Hopf transformation is just W\df — df + W\f = 0. 
For the fc-truncation, we consider the generalization, 



(5.10) 



md k f = + Wl f^ + --- + w k f = o. 

ion under 

-{md k f) 



The invariance of this equation under the evolution (5.8 1 can be shown as follows: 

d_ 



^ - dm d k f + w k d k lf 



(B 3 W k - m®) d k f + Wkd k ®t 



W k d k 



dtj dxi ) 



This implies that the equation W k d k f = is invariant, if / satisfies the linear equations, 



(5.11) 



df frf 



dxi 



for j = 1,2, 



Now we construct the solutions of the KP equation from the linear equations (5.11 1: Let {fj \ j — 
1, . . . , k} be a fundamental set of solutions of (5.10), i.e. 

fj k) + «"i/f _1) + • ■ ■ + ufc-i/j + w k f 3 =0, j = 1, . . . , k . 
Using the Cramer's rule, the function w\ is expressed in terms of the Wronskian with those /j's, i.e., 

d 



Wi 



Ox 



lnWr(A,...,/ fc ). 



We then obtain a solution of the KP equation, 



u = u 2 = 



: w 1 = ^]nWr(f 1 ,...,f k ). 



dx dx 2 

Here the Wronskian is called the r-function of the KP equation (see also [32] )■ Thus the solution 
u in this form is characterized by the kernel of the differential operator, W k d k , and the Wronskian 
structure of the T-function leads to the notion of the Grassmannian as explained below. 

It is now clear that this r-function agrees with r k in (5.2) of the Toda lattice, if we take fj to be 



Si 



.0-1) 



for j=l,2, 



dxi- 1 

One can also show that the square of a k in the L matrix of the symmetric Toda lattice is a solution 
of the KP equation, i.e. 

d 2 



a~ k =u = In r k with r k — Wr(n , t[ , t 



{k-X)s 

i 'i ) ■ 
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Figure 5.2. A 2-soliton solution of the KP equation. The figures show the contour 
plots of the solution u(x,y,t) in the x-y plane with the times for t = —10 (left) 
—5, 0,5 and 10 (right). The r-function is given by Ta(t) in Example 5.1 with 
ti = x,t-2 = y and £3 = t. The parameters (Ai, . . . , A4) are the same as in Figure 



Example 5.2. In Figure 5.2 we show the 2-soliton solution given by ra(t) in Example 5.1 that 
is, u(x,y,t) — 02(^^2,^3) with t\ = x,t% — y and £3 = t. This solution illustrates a resonant 
interaction of two solitons, which was first found in [6]. Notice that there are six distinct regions in 
the ti~t2 plane, four unbounded regions and two bounded regions at the interaction point. Those 
six regions correspond to the dominant exponentials in r 2 function, and those are identified as the 
vertices of the moment polytope of Gr(2,4) (see Section 5.3 and also [IB]). 



5.3. Grassmannian Gr(k,n). There exists a natural identification between the space of the r- 
functions of the Wronskian form Wr(/i, . . . , fk) and the Grassmannian Gr(k,n), the set of k- 
dimensional subspaces in R n . To explain this, we take the following functions as a fundamental 



set of the solutions (5.101, i.e. finite Fourier series solutions of (5.11 1 

n 

(5.12) 



fi 



E 



for 



1. 



, k < n , 



with some constants a^- which define the k x n matrix Arjf . 



,) := (ciij), and the phase functions 9j 



are given by 
(5.13) 



f (t)= E At** for j = l.. 



i=l 



Here A.,- are arbitrary constants which can be identified as the eigenvalues of the L matrix for the 
Toda lattice. We assume A^ being ordered as ( 5.4 1 , i.e. Ai < • • • < A„. 

Since {e Sj : 3 — 1, . . . , n} and {fi : i — 1, . . . , k} are linearly independent sets of functions, it 
follows that Span K {e Sj : 3 = 1, . . . ,n) = M™, and also that Span K {/i : i = 1, . . . , k} defines a k- 
dimensional subspace of R n . We then consider the identifications for the generic set {Ai, . . . , A„}, 



E j := e b 



Vi = l,2, 



and 



Note here that 



M 

:= E a 



K- 

13^3 1 



Vi = l,2, 



,,fc. 



fl j - 1] = {F uej ), j = l,. 

where ej is the j-th vector of the standard basis {ei}™ =1 , and 
a point F G Gr(k,n) is defined by F = Span R {Fi : i = 1, 
represented by the k x n, full rank coefficient matrix A = (oiy) whose rows are the coordinates of 
the basis vectors {Fi : i = 1, . . . , k} with respect to the fixed basis {Ej : j = 1, . . . , n} of K™. Since 



, •) is the scalar product on R n . Then 
,.,k}. Equivalently, F € Gr(k,n) is 



3M 
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A depends on the choice of basis, the matrix representation for the point F £ Gr(k, n) is unique up 
to a left multiplication: A — > HA by any H e GL(k). So, the Grassmannian can be considered as 
the factor space of Mfc Xn (R), the set of all k x n matrices of rank= k, with GL(k), i.e. 

Gr{k,n) = M kxn (R)/GL(k) . 

Note dim Gr(k, n) = k(n — k). The matrix H S GL{k) is chosen so that the coefficient matrix A is 
put in the reduced row-echelon form (RREF) via Gaussian elimination. (Recall that, for an k x n 
matrix in RREF, the leftmost nonvanishing entry in each nonzero row is called a pivot, which is 
normalized to 1. The submatrix of the pivot columns of the full rank matrix A is the k x k identity 
matrix). A convenient paramctrization of Gr(k,n) is given by the Pliicker embedding, 

Gr(k,n) — ► P(A fc M") 
F = Span R {Ji : * = 1, . . . , k] .— > F 1 A — AF k . 

With respect to the basis {E il A ■ ■ ■ A Ei h : i\ < ■ ■ ■ < ik} for A fe K n , 

(5.15) F x A-AFu - £ A(i u ■ ■ ■ , i k ) E it A ■ ■ ■ A E ik , 

1< il <...<ik<n 

where the coefficients A(ii, . . . , ik) of A are the k x k maximal minors of A, called the Pliicker 
coordinates. These are the homogeneous coordinates of the image of F in P(/\ R n ), and they 
satisfy the Pliicker relations, 

fe+i 

E (-ly^Or,^,- . . ,i fe _i) ...,j r ,.. . , jfe+i) = 0, 

r=l 

for any set of numbers {«!,..., ife_i, Ji, . . . , j/'fe+i} S {1,2,..., 2k}. Here j r implies the deletion of the 
index j n . Then Gr(k,n) is also defined as the set of vectors in the form (5.151 with the coordinates 
satisfying the Pliicker relations. Notice that the functions . . . ,i k ) are the Pliicker coordinates, 
that is, they satisfy the Pliicker relations (this is the general fact for the r-functions, see below). 
Let (•, •) be the usual inner product on A k R n , i.e. 

(vx A ■ ■ ■ Av k , w 1 A ■ ■ ■ Aw k ) := det((w i , w 3 ))i< M -< n . 

Then the Wronskian E(ii, . . . ,i k ) is given by 

E(ii, ...,ifc) = (E h A- ■■ A JSi fc ,«[i...,fc]) , 

with W[i v ...fe] — e± A- ■ • Aefc. Note here that the independence of the exponential functions E(i\, . . . , i k ) 
can be shown for the generic choice of the parameters {Ai, . . . , A„}. 
Now let us define the set of r-functions, 

T(k, n) := {t a = Wr(A, . . . , f k ) : (/ 1; . . . , f k ) = (e e \ . . . , e e -)A T ) 

Then comparing the expression of a point on Gr(k,n) and the r-function of the Wronskian deter- 
minant Wr(/i, . . . , f k ), we obtain the following Proposition: 

fe 

Proposition 5.1. Let \(ii, . . . ,i k ) — E ^j%> an< ^ assume that those are all distinct. Then there 

j=1 

exists a bijection, 

T(k,n) Gr(k,n) 

r fe =Wr(A,...,/ fc ) F = F 1 A---AF k 



THE FINITE NON-PERIODIC TODA LATTICE: A GEOMETRIC AND TOPOLOGICAL VIEWPOINT 39 



Proof. Using the Binct-Cauchy theorem, we have 
(5.16) r k = 



E A(ii,...,i k )E(ii,...,i k ) 

l<ii <••• <ifc <n 



where A(i±, . . . , ik) are the Pliicker coordinates given in ( |5.15 1, and 



i k ) = Wr(e e ^,...,e e ^)=l[(X ll - A 4 .) exp 6(i u . . . , i k ) . 

3<l 



Here 6{i\,...,i k ) = E Since ■ ■ ■ , ik)} are distinct, the set {E(i\, . . . , i k )} is linearly 

i=i 

independent as the functions of x. This implies that one can identify the basis {E^ A • ■ ■ A Ei k } of 
A k W l with {E(ii, . . .,i k )} of T(k,n), and we h ave r fc = (Fi A • • • A F k , □ 

We thus identify each r function in the form (5.161 as a point of Gr(k,n), and the solution of the 
KP equation given by T k (t) defines a torus orbit on Gr(k,n), 

[Frit), F k (t)] = [Et(t),..., E k {t)\ A T = Kdiag(e ei « . . . , e 6 ^) A T , 

where K = {hj)i<i,j<n with ki 



Note in particular that 
= C k A E 3 (t) = \*E 3 (t). 



dt k 



where C\ is the companion matrix (3.141. 

As discussed in Section 3.1 let us define the moment map fi : Gr(k, n) — > H^, r k i— ► fi(T k ) with 

(5.17) 



Kn) = 



E |A( il ,..., lfe )e^--^)| 2 (L 4l +--- + i J J 

<•■ ■ <ifc <n 



E |A(*i ifc ) e fl(*x.-.**)|2 

l<ii<"'<^/s<n 



where #(«i, . . . , i k ) = E ^(^* ■ > ano - -^j are the weights of the standard representation of SL(n) 

3 = 1 

(see (3.1 1), and Hjjj is defined by 



7~Lr — Span R < Li, . . . , L n : E A? 



pn— 1 



Then from (5.171, we can see that the image under the moment map of the toric variety generated by 
the solutions of the KP equation with (5.16 1 is a convex polytope whose vertices are the fixed points 
of the orbit. In the representation theory, this polytope is a weight polytop e of the fundamental 
representation of sl(n) on t\ k V with the standard representation V . In Figure 5.3 we illustrate the 
moment polytope of GV(2, 4). The orbit given by a KP solution can be realized as a curve inside of 
the polytope. 



6. The Toda lattice and integral cohomology of real flag manifolds 

Here we briefly explain how one can get the integral cohomology of real flag variety from the 
isospectral variety of the Toda lattice (this is an introduction of the papers [HJUS]): We consider 
the Toda lattice hierarchy (3.3 1 on the real split semi-simple Lie algebra sZ(n,R), and assume X 6 
sl(n, K) to be a generic element in the tridiagonal Hessenberg form, that is, it has all real and distinct 
eigenvalues (see [TBI US] , for the general case associated with real split semisimple Lie algebra) . 
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1100 




0011 



Figure 5.3. The moment polytope of Gr(2, 4). Each vertex corresponds to a fixed 
point of the KP flow given by the r 2 function on R , and it represents the dominant 
exponent in T2, e.g. the weight 1100 = L\ + L2 corresponds to the term e 6l+02 , 
which is dominant as t\ — > —00 for Ai < • • • < A4. 



6.1. The moment polytope and Weyl group action. Let G = SL(n,M,), B be the Borcl 
subgroup of upper triangular matrices of G, and N be the subgoup of lower unipotent matrices of 
G. As in (3.21 1, we consider the companion embedding [TBI [3Tj . 

c A : M A — ► G/B 

X(0) — > n 



(6.1) 



n n 1 mod B 



where n e N is the unique element such that X(0) = noCA^p 1 with the companion matrix (3.141 
[53]. With the factorization (3.4 1 for the hierarchy (3.3), i.e. 

e e(x(o).t) = n ( t)6(t ) with e(x(o), t) = E t x(oy , 

we have 

X(t) = n- 1 (t)X(0)n(t) . 

The Toda lattice hierarchy can be linearized on the flag variety G/B, that is, we have a commutative 
diagram [31] . 

X(Q) - 

A^ri-^t)) j 

X(t) - 

where we have used n ~ 1 n(t) = e s ^ x ^'^ modi? and X(0)rio = uqCa- With this diagram, the 
flows of the Toda lattice hierarchy form a toric variety, denoted by A4a, in the flag manifold G/B, 
that is, 



n 1 mod B 



e e(c A ,t) n -i mod B 



(6.2) 



M A = G^n^B 



ith 



G Ca = { 



= J e 8(C A ,t) 



t e 



In order to characterize this toric variety, let us first consider a diagonal embedding of G/B, 
(6.3) 4> : G/B — ► Gr(l, n) x Gr(2, n) X • • • X Gr(n - 1, n) , 

which defines the Bruhat decomposition, that is, each point in G/B determines a set of points 
(Fi, . . . , V n -x) with Vfe G Gr{k, n) such that 

{0} C Vl C • • • C V n _! C R" , 
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(1,2,3) 




(3,2,1) 



Figure 6.1. The moment polytope for the sZ(3,K.) Toda lattice. The vertices are 
marked by the ordered sets of the eigenvalues, i.e. (i,j,k) — (A*, Aj,Afc). Each 
edge corresponds to a simple reflection rj exchanging the eigenvalues at the jth and 
(j + l)th places. Then the vertices are also marked by the elements of the Weyl 
group, e.g. (1,2,3) = e, (2,1,3) = r\ etc. Each directed edge corresponds to an 
sZ(2,M) Toda flow, and the directed curve in the middle shows a generic flow of 
s/(3,K) Toda lattice. 



where Vk can be expressed in terms of a frame Vk = [Fx, . . . , Fk], where the {Fj : j = 1, . . . , fc} is a 
basis of Vk- As shown in Proposition 5.1 each can be identified with a point on Gr(fc,n) so that 
one may write 

<{> o c A (X(t)) = 0(6^^0 1 B) = (n(t), . . . , 7- w _i(t)) . 



Now we define the moment map /i : G/ B — > 7Yjjj, which is expressed in terms of the map (5.171, 

fi : G/B — ► 7YJj 

(6 ' 4) e^)n^B —> I>(T fc (t)) 

fc=i 

The fixed points of the flows in G/B are the images of the matrices X with all = 0. There are 
n! fixed points, which is the order of the symmetric group S n , and those are the orbit of the Weyl 
group W — S n . As we explained in Section |3.1| we also note that the image under of the toric 

variety M.\ in (6.2) is the convex hull of the vertices corresponding to the fixed points, which are 

j 

the weights L = ikLk with € {0, 1, . . . ,n — 1} and ik ^ ij if k ^ j. The highest weight is 

k=l 

n 

given by = ^ (n — k)Lk, which corresponds to the matrix X(t) as ti — > —oo with the ordering 

fe=i 

Ai < • • • < A„, i.e. diag(X) = diag(Ai, . . . , A„) (see, for example, Figure 3.1 1. 

This can be seen from the dominant exponential term in Tk(t\) for t\ — > — oo, i.e. 

Tk(h) w A(l,...,k)E(l, ...,k) as fi->-oo, 

which implies fi(Tk) = L% + ■ ■ ■ + Lk- Since the vertices are the orbit of the Weyl group W, those 
can be also parametrized by the elements of W = S n . Let r, G If be a simple reflection, that 
is, ri exchanges the entries in positions i and i + 1 on the diagonal of the fixed point matrix. 
Each vertex of the polytope can be also represented by the n-tuple of the eigenvalues, that is, 
(Atu-i( i), A u) -i(2), ■ ■ ■ , A m -i( re )) is the diagonal of the fixed point matrix marked by w <E W. In 
Figure [6T| we illustrate the weight polytope associated with the s/(3, R) Toda lattice where (i, j, k) = 
(A^ Aj,Afc). Each edge given by an arrow corresponds to a simple reflection, and it gives the flow 
associated with sZ(2,R) Toda lattice. 
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6.2. Integral cohomology of G/B. Here we give a brief summary of the cohomology of G/B as 
a background for the next section where we explain how one gets the cohomology of G/B from the 
isospectral variety of the Toda lattice associated with real split semisimple Lie group G. 
Let us first recall the Bruhat decomposition of G/B. 

G/B= |J X w with X w =NwB/B. 

Each Bruhat cell X w is labeled by the element w £ W and codim (X w ) = l(w) where l(w) represents 
the length of w. Then we can define the chain complex, 

C* = C k with C k = % x wi 
fe=o l(w)=k 

where is the longest element of W, and the coboundary operators 5k on X w with l(w) = k is 
given by 

Sk(X w ) = i w '■ w '\ X w , 

where [w : w') is the incidence number associated with X w — ^ X w >. It has been known (see [45| 120]) 
that the incidence number is either or ±2 for the real flag manifold G/B of real split semi-simple 
Lie group G. The Bruhat order defined on W is given as follows: For two elements w and w' in W, 
we have the order 

w < w' iff X w D X w : , 
where the closure X w is called the Schubert manifold. Then the cohomology of G/B can be calcu- 
lated from the incidence graph Qg/b defined as follows: 

Definition 6.1. The incidence graph Qg/b consists of the vertices labeled by w £ W and the edges 
=> defined by 

!(i) w < w' 
(ii) l(w') = l(w) + 1 
(iii) [w : w'] ^ 

The incidence number for each edge is either or ±2 (see [5U]). The integral cohomology is then 
calculated from the graph. 

Example 6.2. In the case of G — SL(3, K), the incidence graph is given by 

e 

[1] [2] 

a- u 

[12] [21] 
[121] 

where the vertices of the hexagon are marked by the elements of the Weyl group denoted as • • • Tj = 
[i ■ ■ ■ j\. Since the nonzero incidence numbers are ±2, the integral cohomology of G/B in this case 
is give by 

H°(G/B,Z) = Z 
H 1 {G/B 1 %) = 



H 2 (G/B,Z) = Z 2 ®I 
H 3 (G/B,Z) = Z 



'2 



In Figure 6.2 we show the incidence graph for G/B with G — SL(A, M), from which one can 
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[2] [1] P] 

[23] [21] [12] [13] [32] 

[213] [121] [232] [321] [123] [132] 

[1231] [2132] [2321] [123 2] [1321] 

[12312] [23212] [12321] 

w ,.=[123121] 



FIGURE 6.2. The incidence graph Gg/b for the real flag manifold SL(A,R)/B. The 
Bruhat cells X w — NwB/B are denoted by [ij...k] for w — r^j . . .r^. w* — 
[123121] is the longest element of W = S4. The incidence numbers associated with 
the edges are ±2 (see also Example 8.1 in [20 ). 



compute the integral cohomology [2D]. We then obtain 



H°(G/B,Z 
H l {G/B,Z 
H 2 (G/B,Z 
H 3 (G/B,Z 
H 4 (G/B,Z 
H 5 (G/B,Z 
H 6 (G/B,Z 



Z 


Z 2 © Z 2 ffi z 2 

z e z © z 2 e z 2 
z 2 ez 2 
z 2 © z 2 ® z 2 
z 



The incidence graph for the general case of real split semisimple G can be found in 20J . Then the 
integral cohomology of G/B can be computed from the incidence graph with the incidence numbers 
[w : w') being or ±2. 

For the rational cohomology, we have 

(6.5) H* (G/B, Q) = if* (K/T, Q) = H* (K, Q) , 

where i"T is the maximal compact subgroup of G, and T is the maximal torus of G. e.g. for 
G = SX(n, E), K = SO(n), and T = diag(±l, . . . , ±1) (see Proposition 6.3 in Q2]). It is well known 
that the cohomology ring H* (K, Q) of compact connected group K of rank / is given by the exterior 
product algebra, 

H* (K, Q) =/\ (tj{a; mi , x m2 ,...,x mi }, 

where {x mi 



, x mi } are the generators of the exterior product representation with deg(x TOi ) 



, I and 



mi 



■mi = dim(K) (see for example [14 ). In the case of K = SO(n), 



(odd) for i = 1, 
we have, 

(a) for n = 2m + 1, 

H*(SO(2m + 1),Q) = f\q(x 3 ,x 7 , . . .,x 4m _i) 

(b) for n = 2m, 

H*(SO(2m),Q) = f\®(x 3 ,x 7 , . . . , x 4m _ 5 , y 2m -i) ■ 

Note here that the generators include the additional y 2m _i. For example, i?*(50(4),Q) is 
generated by two elements a; 3 , y 3 of the same degree, V — Qx 3 + Qy 3 and h 2 V = Qx 3 A y 3 . 
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We also note that the number of points on the finite Chevalley group K(¥ q ) of the compact 
connected group K is given by certain polynomial of q. Here ¥ q is a finite field with q elements. 
Although this polynomial can be computed by using the Lefschetz fixed point theorem for the 
Frobenius map <f> : K(¥ q ) — > K(¥ q ), x i— > x q , we here give an elementary calculation to find those 
polynomials for K — SO(n) (see also |18j). As we show in the next section, those polynomials are 
also related to the indefinite Toda lattice. 

Let us first assume that q is a power of a prime number p ^ 2, such that in ¥ q the polynomial 
x 2 + 1 is not irreducible, i.e. \/— T € F 9 . Then we have the following results for |£ n (F 9 )|, the number 
of ¥ q points on S n : 



\S n (¥ q )\ = 



Lemma 6.1. We have 

q m - 1 (q m - 1) if n= 2m- 1, 
q m (q m + 1) if n= 2m. 

Proof. Let us first consider the case n = 1, i.e. 

S 1 (¥ q ) = {(x,y)e¥ 2 q :x 2 + y 2 = l}. 

Then using the formulae for the stereographic projection; x = M 2 " x , y = ^ 2 ~* with y ^ 1 and 
{« £ F, : it 2 + 1 ^ 0}. Since \/— 1 G ¥ q , we have 2 points in {u 2 + 1 = 0}. Counting the point 
(0, 1), the north pole, we have |S ,1 (F g )| = q — 2 + 1 = q — 1. Now consider the case n = 2, we have 

* = S^+T.W = i^+T.* = S^f+T with ^ 1 and {(m,u 2 ) £F>? + ^ + M0}. This gives 
g 2 — (<? — 1) points (note (g — 1) is the number of points in u\ + + 1 =0). We now add the points of 
the north pole (x, y, 1) with x 2 + y 2 = 0. This gives 2(q — 1) + 1, where 2(q — 1) for x = ±\/—ly ^ 
and 1 for (0,0,1). Then we have |5 2 (F,)| = q 2 - (q- 1) + 2(q- 1) + 1 = q(q+l). Using the induction, 
one can show that the number of points in the north pole is given by |{(xi, . . . ,X2 m -i) & F 2m_1 : 
3!? + - ' -+^lm-i = 0}| = q 2m ~ 2 and \{x u . . .,x 2m ) € F 2m : x\ + - ■ - + x 2 2m = 0}| = q>m-i +q m _ q m-i _ 
Then one can obtain the above formulae for |5"(F q )|. □ 

We can now find the number of ¥ q points of finite Chevalley group SO(n, ¥ q ): First recall that 
SO(n + l)/SO(n) S n . Then we obtain 



\SO(n,¥ q )\ = Y[\S n ' k (¥ 



which leads to the results [T^] : 

(a) For n = 2m, 

\SO(2m,¥ g )\ = 2q m{m - 1 \q 2 - l)(g 4 - 1) ■ • • (q 2m - 2 - l)(q m - 1) . 

(b) For n = 2m + 1 , 

\SO(2m + l,¥ q )\ = 2q m \q 2 - 1) (g 4 - 1) • • • (q 2m - 1) . 

In general, the number of ¥ q points on the compact group K can be expressed by (see e.g. |14p 

i 

(6.6) \K(¥ q )\=q r p(q) with p(q) = JJ («* - 1) , 



where di's are degree of basic Weyl group invariant polynomials for K given by di — (rrij + 1) /2, and 
r = dhn(K) — deg(p(g)). In the next section, we show that those polynomials can be reproduced by 
counting the blow-ups in the solution of the indefinite Toda lattice (see [TB] for the general case) . 
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6.3. Blow-ups of the indefinite Toda lattice on G and the cohomology of G/B. Now we 
explain how we get the cohomology of G/B from the moment polytope of the indefinite Toda lattice 
of Section \2A\ 

First we note that the r-functions can change their signs if some (but not all) of Sj's are negative. 
This can be seen from (2.40 1, and implies that the solution blows up for some time t\ = t\, (see also 



(2.281 and ( 2.29 1 ) . The explicit form of the r-functions can be obtained from (2.271, and they are 



expressed by (see (5.16) and also Proposition 3.1 in [5T]1. 



(6.7) 



T k (t) 



E 



i<ji<--<jfc<™ 



s ik K (h, ■ ■ -,3k) exp 



E A A t 



where K(ji, . . . ,jk) are positive and given by 

K(j 1 ,...,j k )=( V >°(\ h )---<p°(\ jk )y : . : .,.). 
As the simplest case, let us consider the sl(2,R) Toda lattice: We have one r-function, 



n(t) = s lPl e Xlt + s 2 p 2 e X2t . 

If sis 2 < 0, Ti(t) has zero at a time t = Aq 1 Ai ln(^), that is, we have a blow-up in the solution. 
The image of the moment map /i(ri) is given by a line segment whose end points correspond to 
the weights L\ and L 2 — —L\. Although the dynamics are so different in the cases s\s 2 > and 
si-s 2 < 0, the moment polytope (a line segment) is independent of the signs of the s,'s. Notice that 
si-s 2 — sgn(gi), and in general, if sgn(gfe) < for some k, then the solution blows up sometime in M. 

In order to find the general pattern of the sign changes in (31 (t), . . . , <?„_i(i)) of the matrix X 
in (2.221, we first recall that the isospectral variety is characterized by the moment polytope M. t 
whose vertices are given by the orbit of Weyl group action. Here the set of signs e = (ei, . . . , e n -i) 
is defined by the signs of gi for t — ► —00. From the ordering Ai < • • • < A n , we first see that 
s k K(l,. . . ,k)exp((X 1 + ■ 



r k (t) 



si 



Afc)i). Then from the definition of gk(t) in (2.281, i.e 



9k = T k -iT k+ i/T k , the sign of g k {t) for t -> -00 is given by 

e k = sgn(5 fe ) = s k s k+ i for k = 1, 



1 



Then from the moment map (5.171, one notes that the moment polytope given as the image of the 



moment map in (6.4 1 is independent of the sign set e. However the dynamics of the Toda 

lattice with a different e is quite different, and the solution with at least one e k < has a blow-up 
at some fef. 

We now consider each edge of the polytope which corresponds to an s/(2,]R) indefinite Toda 
lattice, that is, where gj 7^ for only one j. This edge can be also expressed by a simple reflection 
Tj G W. Since the simple reflection rj exchanges Sj and Sj+i, we have an action of Tj on all the 
signs efc, Tj : e k -> e' k , 




if 
if 
if 



j = k-l 
j = k + 1 
j = k, or |j 



- k\ > 1 

This formula can be extended 



which can be also shown directly from the form of r k (t) in ( |6.7 
to the indefinite Toda lattice on any real split semisimple Lie algebras, and we have (see (4.4 1 and 
Proposition 3.16 in [15]): 
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Proposition 6.1. Let e 3 = sgn(<7j) for j = 1, . . . , n — 1. Then the Weyl group action on the signs 
is given by 

—Ckj 

where (Cij)i<ij< rl _i is the Cartan matrix o/ sZ(ti, IR). 

With this VF-action on the signs e = (ei, . . . , e„_i) with = sgn(gfc) at each vertex of the 
polytope, we now define the relation between the vertices labeled by w and w' = wri as follows: 
Notice that if ti — +, then (ei, • • • , e„_i) remains the same under reaction. Then we write 

w =>• w' with w = wri . 

Now the following definition gives the number of blow-ups in the Toda orbit from the top vertex e 
to the vertex labeled by w G W: Choose a reduced expression w = rj t ■ ■ ■ rj k . Then consider the 
sequence of the signs at the orbit given by w-action, 

We then define the function 77(11;, e) as the number of — > which are not of the form =>. The number 
77(11;* , e) for the longest element gives the total number of blow-ups along the Toda flow in the 
polytope of A4 e . Whenever e = (— , . . . , — ), we just denote 77(11), e) — rj(w). This number 77(11;, e) 
does not depend on the choice of the reduced expression of w (see Corollary 5.2 in [IB])- Hence 
the number of blow-up points along the trajectories in the edges of the polytope is independent of 
the trajectory parametrized by the reduced expression. In Figure |2.3[ we illustrate the numbers 

i](w,e) for the sZ(3,R) indefinite Toda lattice. For example, on M. , we have 77(e) = 0,77(7-1) = 

vi 1 *?) = v{ r i r 2) = 17(^2^1 ) = 1 and rj{rir2r\) = 2, i.e the total number of blow-ups is 2. We also 
illustrate this for the sZ(4,R) Toda lattice in Figure 6.3 Along the path shown in this Figure, we 
have 77(e) = 0,?7([2]) = jy([21]) = t?([213]) = 1,t/([2132]) = 2, jj([21323]) = 3 and 77(11;,) = 4, where 
[ij ■■■k} = nrj ■ ■ ■ r k , and note [21323] = [12312]. 

In general, the total number of blow-ups 77(111*, e) depends only the initial signs e — (ei, . . . , e„_i) 
with ti = sgn(<7i(r;)) for t — > —00, which is given by ej = SjSj+i. Then in the case of sl(n,M) 
indefinite Toda lattice, the number 77(11^,6) = 771(71 — m) where m is the total number of negative 
s^s (Proposition 3.3 in |51j). In particular, the maximum number of blow-ups occurs the case with 
e =(—,..., —), and it is given by [(71 + l)/2](n — [(71 + l)/2]). Those numbers 77(11;*, e) are related 
to the polynomials appearing in ¥ q points on certain compact groups defined in (6.6 1. 

We now introduce polynomials in terms of the numbers 77(11;, e), which play a key role for counting 
the number of blow-ups and give a surprising connection to the rational cohomology of the maximal 
compact subgroup SO(n) (Definition 3.1 in [T8]). 



Definition 6.3. We define a monic polynomial associated to the polytope A4 e , 

K<Z,e) = (-l)^*> E (-l) l ^q^\ 

where l(w) indicates the length of in. Notice that the degree of p(q, e), denoted by deg(p(q,e)), is 
the total number of blow-ups, i.e. rj(w*, e) = deg(p(g, e)). For the case e = (— , ...,—), we simply 
denote it by p(q) . 

Example 6.4. In the case of the sZ(2,R) Toda lattice, 

(a) for e = (+), we have e => si which gives p(q,+) = 0, 

(b) for e = (— ), we have a blow-up between e and si, hence p(q, — ) = q — 1. 

Recall from the previous section that the polynomial p(q) = p(q, — ) appears in 15*0(2, ¥ q ) \ = q — 1. 
In the case of the sZ(3,R) Toda lattice, from Figure [273) 
(a) for all the cases of e = (ei, e 2 ) except (—,—), we have p(q, e) = 0. 



THE FINITE NON-PERIODIC TODA LATTICE: A GEOMETRIC AND TOPOLOGICAL VIEWPOINT 47 




FIGURE 6.3. The moment polytope M for the s/(4,R) indefinite Toda lattice. 

The divisors defined by the set of zero points for the r-functions are shown by the 
dotted curve for {ti — 0}, by the light color one for {r 2 = 0} and the dark one 
for {T3 = 0}. The double circles indicate the divisors with {r^ = 0} n {tj = 0} 
which are all connected at the center of the polytope corresponding to the point 
with {rj = r 2 = r 3 = 0}. The numbers in the polytope indicate the number of 
blow-ups along the flow. An example of a path from the top vertex e to the bottom 
vertex w*, the longest element of S4, is shown by directed edges. 



(b) for e — (— , — ), we have p(q) — q 2 — 1. 
Note again that the polynomial p(q) appears in |50(3,F 9 )| = q(q 2 — 1). 
In the case of sZ(4, R), we have, from Figure [6~3) 

(a) for all e= (ei,e 2 ,e 3 ) except (-,-,-), p(q,e) = 0. 

(b) for e = (-, - -), p(q) = q i -2q 2 + l = (q 2 - l) 2 . 

Again note that \SO(4,¥ q )\ = q 2 (q 2 - l) 2 . 

Casian and Kodama then prove that the polynomial p{q) for M. e with e = (—,...,—) in Definition 



6.3 



agrees with the polynomial p(q) in |lf(F g )| in (6.6) where K is the maximal compact subgroup 
of real split semisimple Lie group G for the Toda lattice (Theorem 6.5 in [T5]V 

Thus the polynomial p(q) contains all the information on the ¥ q points on the compact subgroup 



K of G, which is also related to the rational cohomology, i.e. H*(K,Q) — H*(G/B,Q) (see (6.5)) 



Now recall that the integral cohomology of the real flag variety G/B is obtained by the incidence 
graph Q G i b in Definition|6.1| In |18j , Casian and Kodama show that the graph Q G j B can be obtained 
from the blow-ups of the Toda flow. They define a graph Q t associated to the blow-ups as follows: 

Definition 6.5. The graph Q e consists of vertices labeled by the elements of the Weyl group W and 
oriented edges =>. The edges are defined as follows: 

(a) w\ < W2 (Bruhat order) 

(b) l( Wl ) = l(w 2 ) + 1 

(c) rj(w 1 ,e) = r)(w 2 ,e) 

(d) w^e = w^ x t 
When e = (— , . . . , — ), we simply denote Q = Q e . 

Then they prove that Q e with e = ( — , ...,—) is equivalent to Gg/b (Theorem 3.5 in |18| which is 
the main theorem in the paper) . For example, the graph Q associated with Figure |6.3| agrees with 



W\ =>■ w 2 iff 
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the incidence graph Qg/b given in Figure 6.2 The proof of the equivalence Qg/b — Q contains 



several technical steps, which are beyond the scope of this review. 
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